TOWARD REGIONAL CHARACTERIZATIONS 
OF THE OCEANIC INTERNAL WAVEFIELD 
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Many major oceanographic internal wave observa- 
tional programs of the last four decades are reanalyzed 
in order to characterize variability of the deep ocean 
internal wavefield. The observations are discussed in 
the context of the universal spectral model proposed by 
Garrett and Munk. The Garrett and Munk model is a 
good description of wintertime conditions at Site-D on 
the continental rise north of the Gulf Stream. Elsewhere 
and at other times, significant deviations in terms of 
amplitude, separability of the 2-D vertical wavenumber 
- frequency spectrum, and departure from the model's 
functional form are reported. 

Specifically, the Garrett and Munk model overesti- 
mates annual average frequency domain spectral levels 
both at Site-D and in general. The bias at Site-D is 
associated with the Garrett and Munk model being a 
fit to winter time data from Site-D and the presence of 
an annual cycle in high frequency energy in the western 
subtropical North Atlantic having a maximum in win- 
ter. The wave spectrum is generally non-separable, with 
near-inertial waves typically having greater bandwidth 
(occupying smaller vertical scales) than continuum fre- 
quency waves. Separability is a better approximation 
for more energetic states, such as wintertime conditions 
at Sitc-D. Subtle geographic differences from the high 
frequency and high vertical wavenumber power laws of 
the Garrett and Munk spectrum are apparent. Such de- 
viations tend to co-vary: whiter frequency spectra are 
partnered with redder vertical wavenumber spectra. 

We review a general theoretical framework of sta- 
tistical radiative balance equations and interpret the 
observed variability in terms of the interplay between 
generation, propagation and nonlinearity. First: Non- 



linearity is a fundamental organizing principle in this 
work. The observed power laws lie close to the In- 
duced Diffusion stationary states of the resonant kinetic 
equation describing the lowest order nonlinear transfers. 
Second: Eddy variability and by implication wave-mean 
interactions are also an organizing principle. Observa- 
tions from regions of low eddy variability tend to be 
outliers in terms of their parametric spectral represen- 
tation; other data tend to cluster in two regions of pa- 
rameter space. More tentatively, the seasonal cycle of 
high frequency energy is in phase with the near-inertial 
seasonal cycle in regions of significant eddy variability. 
In regions of low eddy variability, the seasonal cycle in 
high frequency energy lags that of near-inertial energy. 
The Induced Diffusion stationary states are approx- 
imate analytic solutions to the resonant kinetic equa- 
tion and the Garrett and Munk spectrum represents 
one such analytic solution. We present numerical solu- 
tions of the resonant kinetic equation, though, that are 
inconsistent with the Garrett and Munk model repre- 
senting a stationary state, either alone or in combina- 
tion with other physical mechanisms. We believe this 
to be the case for other regional characterizations as 
well. We argue that nonstationarity of the numerical 
solutions is related to local transfers in the horizontal 
wavenumber domain whereas the analytic Induced Dif- 
fusion solutions consider only nonlocal transfers in the 
vertical wavenumber domain. Consequences for under- 
standing the pathways by which energy is transfered 
from sources to sinks are considered. Further progress 
likely requires self-consistent solutions to a broadened 
kinetic equation. 



1. INTRODUCTION nism for transferring energy and momentum across 

large distances and across different scales. In par- 
Internal waves are ubiquitous features oi geo- 
physical fluids, contain a significant fraction of the ticular, the drag associated with internal wave 

total variance and constitute an effective mecha- breaking needs to be included in order to obtain ac- 



Polzin and Lvov 



pages 1-57 



• 1 • 



2 • Polzin and Lvov: REGIONAL CHARACTERIZATIONS 



curate simulations of the atmospheric Jet Stream 
[Fritts et al., 2003] and it has been argued that the 
ocean's Meridional Overturning Circulation [Wun- 
sch and Ferrari, 2004] is forced by the diffusion of 
mass [Ledwell et al., 2000] associated with internal 
wave breaking [Polzin et al., 1997] rather than by 
the convective production of cold, dense water at 
high latitudes. Both circulations represent funda- 
mental pieces of the Earth's climate system. 

The oceanic internal wavefield is "complex". It 
spans a vast range of scales: in the time domain 
periods range from about one day for inertial waves 
to a few minutes for buoyancy frequency waves in 
the pycnocline. Horizontal wavelengths range from 
a few meters to hundreds of kilometers. The inter- 
nal wave field's complexity arises not just from its 
extended range of scales, but also from a multiplic- 
ity of possible sources at the surface (atmospheric 
forcing associated with wind stress and buoyancy 
fluxes) and bottom (tides and mean currents im- 
pinging upon nonuniform topography); the poten- 
tial for scale transformations (reflection and scat- 
tering from topography) and dissipation at the bot- 
tom boundary; and interior transfers of energy and 
momentum through wave-mean flow interactions. 
A significant sink is wave breaking, which is be- 
lieved to be product of nonlinearity resulting in 
downscale transfers of energy. 

Despite, or perhaps because of, this complex- 
ity, a universal character is usually ascribed to 
the oceanic internal wave energy spectrum. In 
a classic work, Garrett and Munk [1972] demon- 
strated how observations from various sensor types 
could be synthesized into a combined wavenumber- 
frequency spectrum, now called the Garrett-and- 
Munk (GM) spectrum of internal waves. Con- 
sistent only with linear internal wave kinematics, 
the GM spectrum was developed as an empirical 
curve fit to available data. Deviations have been 
noted near boundaries [Wunsch and Webb, 1979; 
Polzin, 2004b] and at the equator [Eriksen, 1985a]. 
However, our notion of what constitutes the back- 
ground wavefield has remained static since the last 
significant model revision provided by Cairns and 
Wilhams [1976]. 

The authors' sentiments are captured by the 
first paragraph of Briscoe [1975a] 's Survey sec- 
tion: "The single paper motivating the most com- 
ments, experiments, and disquiet in a lot of read- 
ers was Garrett and Munk [1972], now called fa- 
miliarly GM72. (The updated version is Garrett 
and Munk [1975], cahed GM75.) The paper is a 
virtuoso orchestration of synthesis, approximation, 
boldness, normalization, and implication. Start- 
ing from the observation that the wavenumber- 
frequency (w — «;) spectrum of a multimodal dis- 
persive system is not uniquely determined by pure 
wavenumber and pure frequency spectra, the pa- 



per goes onto construct from linear dynamics and 
WKBJ methods anoj — K spectrum that is astonish- 
ingly consistent with many kinds of data from all 
kinds of instruments. Wunsch [1975] iconoclasti- 
cally discusses the pitfalls of such an approach and 
correctly emphasizes the fact that the dynamically 
important trait is the inconsistency of various data 
sets with some base state model." With immense 
inspiration provided by Garrett and Munk, we take 
up, after 35 years, the difficult task posed by Wun- 
sch. We document variability of the oceanic inter- 
nal wavefield and then review the field to assess the 
dynamical underpinnings of that variability. While 
heavily referenced, this work does not attempt to 
be a literature review. 

In previous work [Lvov et al., 2004] we noted 
that deviations of frequency/vertical wavenumber 
power laws from the universal model tend to co- 
vary: whiter frequency spectra are partnered with 
redder vertical wavenumber spectra. Here we pro- 
vide a more systematic study of such deviations. 
We define spatial and temporal patterns of spectral 
parameters (amplitude, power laws, bandwidth 
and separability) in a review of many observational 
programs of the last four decades. Having doc- 
umented variability in the background wavefield, 
we interpret the spatial patterns in terms of gen- 
eration, nonlinearity, wave-mean interactions and 
dissipation. 

Nonlinearity plays a key role. In a compan- 
ion manuscript [Lvov et al., 2010] we discuss the 
application of wave turbulence formalisms to this 
problem and rigorously demonstrate the validity 
of conventional wisdom regarding a truncated an- 
alytic approximation to nonlinear spectral trans- 
fers described in McComas and Miiller [1981a, b]. 
The observed power law combinations in frequency 
and vertical wavenumber documented here are 
largely consistent with the stationary states associ- 
ated with the truncated analytic approximations. 
While this appears to be a stunning affirmation of 
this approximation to the wave turbulence trans- 
fer integral, it comes with caveats, disclaimers and 
possible contradictions. First, the theory is not for- 
mally valid at high wavenumber for realistic wave 
amplitudes. Second, the identification of GM76 be- 
ing a stationary state is only true after integration 
over the frequency domain and invocation of dissi- 
pation. Finally, numerical evaluations of the non- 
linear transfers are inconsistent with the high fre- 
quency, high vertical wavenumber domain of GM76 
being a stationary state. This tension also applies 
to other stationary states identified by the trun- 
cated analytic theory having power law combina- 
tions associated with energy transfers from larger 
to smaller vertical scales and from higher to lower 
frequencies. Such stationary states require an en- 
ergy source at high frequencies. But our review 
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of the remaining generation and spectral transport 
mechanisms does not yet reveal what this mecha- 
nism might be. 

The article is organized as follows. In Section 2 
we summarize the development of the GM model 
through its various incarnations. We then re- 
view and re-analyze data obtained by major inter- 
nal wave observational programs of the last four 
decades in Section 3. These data clearly demon- 
strate that deviations from the GM spectrum are 
not random, but rather form a distinct pattern 
which we summarize in Section 4. We discuss the 
theoretical framework which may be used to in- 
terpret these observations in Section 5 and con- 
clude in Section 6. An appendix contains addi- 
tional useful information relating to the interpre- 
tation of oceanographic data. 



2. BACKGROUND 

2.1. Kinematic Structure 
2.1.1. Dispersion Relation 

Internal waves arise through the restoring force 

of gravity in a stably stratified fluid. In a rotating 

system they comprise a 'fast' mode of oscillation 

with frequencies between inertial / and buoyancy 

A^ frequencies. In terms of an Eulerian description 

with velocity vector u = {u,v,w), the linearized 

equations of motions are 



ut- fv = - 


-T^x 


vt + fu = - 


-% 


Wt = -T^z 


-h 


Pt + WPz 


= 



Wx + -Viy + Wz 







in which kinematic pressure vr = "pj Po and buoy- 
ancy h = —gp/po are perturbations from a hydro- 
statically balanced equilibrium state, Ifz = — &, 9 is 
gravity and the fluid is assumed to be incompress- 
ible (V • u = 0). We use r = {x,y,z) as a posi- 
tion vector. An over-bar {(/)) indicates a mean on 
time scales much longer than a wave period. The 
constant po is the mean density and arises from 
assuming density to be constant in computing mo- 
mentum changes associated with accelerations, but 
taking density changes into account when they give 
rise to buoyancy forces (the Boussinesq approxima- 
tion) . The buoyancy gradient relates to oscillations 

of frequency N = Jb^ = \/ —gpl Po corresponding 
to periods of minutes in the upper ocean to hours 
in the abyss. This work is concerned almost ex- 
clusively with horizontal scales small enough that 
the Earth's rate of rotation about the local ver- 
tical, / = 2r2sin(latitude), can be considered as 
constant (the f-plane approximation). The rate of 
rotation of the Earth about its axis is 0. The Cori- 



olis frequency / thus corresponds to periods of 1/2 
day at the poles to several days in the tropics and 
vanishes at the equator. The foundation for inter- 
preting observations lies in linear polarization rela- 
tions (Section 2.1.2) obtained from these equations 
of motion. 

The equations of motion can be manipulated to 
obtain a wave equation in vertical velocity (e.g.. 
Gin [1982]): 

(dl + f)dlw + {N\z) + d}\{dl + dl)w = (1) 



for arbitrary stratification profile A^ 



Assum- 



ing the background stratification is constant, plane 
wave solutions w = Woe^'-'^'^~'^^' with three dimen- 
sional wavevector p = {k,l,m) having vertical 
component m and horizontal magnitude fc/^ =| k | 
return 



CTr 



2f2 + j^2q 



m 



m? + kf^ 



(2) 



which can be rewritten as 

.2 f2 



a 



f 



A2 



a^ 



A-2 



Many of the important properties of internal 
waves are immediately available: 

• Incompressibility (V-u = 0) and a plane wave 
solution imply particle velocities are normal to the 
wavevector, p • u = 0. 

• Freely propagating waves, for which p is real, 
require / <| c |< A. 

• Wave frequency solely determines the angle of 
the wavevector p relative to the vertical and, since 
p-u = 0, the ratio of horizontal to vertical velocity. 

• The group velocity Cg = Vpcr (= [ak/k\, 
al/kl, -a/m] if f ^ a^ ^ N^) has the 'odd' 
property that crests and troughs (i.e., lines of con- 
stant phase) propagate in the same direction as 
a wave packet (group) in the horizontal but the 
signs of phase and group velocity differ in the verti- 
cal coordinate. Upward phase propagation implies 
downward energy propagation. 

2.1.2. Polarization Relations 

A suite of diagnostic relations can be obtained by 
assuming constant stratification rate, plane wave 
solutions in which all variables are proportional to 
gi(pr-o-t)g^j^^ expressing the dependent variables in 
terms of wave amplitude a: 
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The prefactor in these polarization relations is such 
that the wave amplitude a is normalized to repre- 
sent the total energy: 

Ek + E^ + Ep = aa*. (8) 

The total energy is the sum of horizontal kinetic 
-Efc, vertical kinetic Ey and potential Ep = i^N'^rj^ 
energy in which r/ represents the vertical displace- 
ment. 

Two important diagnostics utilized here are: 
First, assuming constant stratification, the ratio of 
horizontal kinetic to potential energy for a single 
internal wave is: 



^2 + /2 iV2 



Ar2 



(9) 



Ek uu* + vv 

~Ep " bb*N~^ ^ ^2 - /■ 

Second, the horizontal velocity trace for high fre- 
quency internal waves is rectilinearly polarized. 
Rotation alters this so that the velocity trace at in- 
ertial frequency inscribes a circle with time. With 
depth, the near inertial velocity vector traces out 
an elliptical helix in which the ratio of velocity vari- 
ance in the clockwise (Ecw) and counter-clockwise 
(Eccw) rotating components is given by: 

Ecw ^ {u-iv){u- iv)* ^ (cr + /)2 
Eccw (u + iv) (u + iv)* (cr - ly 

with the sign convention that positive frequency 
f < a < N implies upward phase propagation (and 
hence downward energy propagation) in the north- 
ern hemisphere. 
2.1.3. Weak Spatial Inhomogeneities 

If the stratification profile varies much more 
slowly than the wave phase, a WKB approxima- 
tion for vertically propagating waves provides the 
approximate solution to (1): 

, ^'^^ fNiz)dz 
u;ociV(z)-^/2gV^2_^2 7 ^^^^ 

and so the effects of a variable buoyancy profile 
can be accounted for by stretching the depth co- 
ordinate by A^. This amounts to scaling the ver- 
tical wavenumber by N/No, in which Nq is a ref- 
erence stratification. The value No = 3 cph is of- 
ten used. The buoyancy scaling of other depen- 
dent variables follows from the polarization rela- 
tions (3)- (7). The use of the WKB approximation 
requires a'^ <C N'^- If this relation is not satisfied, 
solutions can be found by treating (1) as an eigen- 
value problem with appropriate boundary condi- 
tions. If (t2 <c A^2 but the wave-phase is not slowly 
varying, the boundary conditions are that w = at 



the top {z = H) and bottom {z = 0), which then 
implies the horizontal velocities [u{z),v{z)] are pro- 
portional to: 

N{zf/'^cos{n^ r N{z)dz/ f N{z')dz')il2) 
Jo Jo 

for integer values of n. 

Similarly, one can investigate interactions of the 
internal wavefield with a mean velocity field by as- 
suming the internal wavefield has small amplitude 
and small spatial scales relative to a time inde- 
pendent geostrophically balanced fiow, Section 5.2. 
This problem is amenable to WKB and ray-tracing 
techniques. In this slowly varying limit, the linear 
polarization relations reported here are modified 
by replacing the Eulerian frequency a with the in- 
trinsic frequency lo = a — p u. 
2.1.4. Slow Oscillations 

Internal waves are to be distinguished from a 
'slow' mode consisting of, at lowest order, the 
steady geostrophic balance. 

1 
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Higher order contributions are time dependence, 
nonlinearity and the effects of a variable rate of 
planetary rotation (/). Substitution of a plane 
wave solution into the geostrophic balance returns 
primarily the information that small aspect ratio 
fluctuations {kh/m <^ f/N) have E^/Ep ratios of 
less than one. 

There are two main paradigms for considering 
these slow motions vis-a-vis the fast internal wave- 
field. The first, already alluded to, is the limit that 
geostrophic velocities are larger and have signifi- 
cantly larger spatial scales. The second paradigm 
comes about as the two share similar spatial scales. 
In this instance the slow motions have typically 
smaller amplitudes and hence appear as a 'con- 
tamination' to the wavefield. Beyond the issue of 
aspect ratio, the consistency relations for this non- 
propagating mode have little utility and this 'con- 
tamination field' has the character of uncorrelated 
noise, which is the content of several diagnostic 
studies [D'Asaro and Morehead, 1991; Kunze and 
Sanford, 1993]. The characterization of this slow 
mode as 'noise' should not surprise. The dynamics 
in this limit are those of stratified rotating turbu- 
lence in which motions at different isopycnals are 
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uncoupled at scales larger than those over which 
shear instabilities can develop [Riley and Lelong, 
2000]. 

While diagnostic models of the slow mode usu- 
ally assume an incoherent field [M tiller et al., 1978], 
the interaction of the internal wavefield with quasi- 
permanent density finestructure can be exploited 
by assuming the wavefield buoyancy scales on the 
density finestructure [Polzin et al., 2003]. This 
buoyancy scaling results in a correlation between 
internal wave shear and quasi-permanent density 
gradients when such a correlation is inconsistent 
with fields composed of either pure waves or pure 
finestructure. This enables an estimate of the 
quasi-permanent spectrum, Polzin et al. [2003] and 
Section 3.3.4. 

One finestructure consistency relation has great 
importance. The plane wave solutions of inter- 
nal waves have no Ertel potential vorticity signa- 
ture, whereas the non-waves have non-zero Ertel 
potential vorticity [Miiller et al., 1986]. The dis- 
tinction is important for issues of isopycnal disper- 
sion [Polzin and Ferrari, 2004] and for dynamics 
[Miiller, 1995]. 

Reliance upon the plane wave solutions, how- 
ever, has its limitations. While internal waves 
and slow modes can be distinguished in terms of 
Ertel potential vorticity content, this distinction 
disappears when one considers a slowly varying 
wavepacket with spatial structure in both hori- 
zontal directions, i.e., solutions of the form 

See Biihler and Mclntyre [2005] and Polzin [2008a] 
for further details. 

2.2. Parametric Spectral Representations 

Much of the original Garrett and Munk [1972] 
paper is a demonstration of how observations from 
various sensor types could be synthesized into a 
combined wavenumber-frequency spectrum consis- 
tent with linear internal wave kinematics. The 
toolbox contains: 

• a linear internal wave dispersion relation (2) 
to transfer from one domain to another, e.g. 

E(kfi)dkfi = E{m)dm 

in which dkh/dm follows from (2), 

• diagnostic relations (Section 2.1.2) to test 
consistency with linear wave kinematics, and 

• the choice of a vertical spectral rather than 
modal representation implies the invocation of an 
'equivalent continuum' in which the waves are as- 
sumed to be vertically propagating. Boundary con- 
ditions and turning points are neglected. See Gar- 
rett and Munk [1972] for a discussion of modal 
properties. 



Assuming that the buoyancy profile varies much 
more slowly than the wave phase implies the to- 
tal energy E, obtained as an integral of the energy 
density over the spectral domain. 



£'(k, cr)(ik(icr = / E{m,a)dmda (13) 



varies as 



E = Er^ 



N 
'Wo' 



(14) 



i.e., E/N is an adiabatic invariant. The energy 
spectrum is denoted as E with following arguments 
[e.g., E{a)]. The variable E without following ar- 
guments represents total energy.) The factors E^ 
and No represent reference values for the total en- 
ergy of the internal wavefield and stratification. 

Garrett and Munk proposed that the spectral 
energy density can be represented as a separa- 
ble function, i.e., the product of a function A of 
vertical wavenumber only and function B of wave 
frequency only: 



E{a, m) 



EA{—)B{a) 



(15) 



where m* is a fixed reference number. This refer- 
ence wavenumber is conveniently expressed as 



W-* = TTJ^/b, 



(16) 



in which the variable j represents the mode num- 
ber of an ocean with an exponential buoyancy 
frequency profile having a scale height of 6 [6 = 
1300 m in the GM model]. Separability in verti- 
cal wavenumber and frequency was invoked as the 
simplest representation not inconsistent with avail- 
able data. The spectral amplitude was regarded 
as being independent of spatial coordinate apart 
from the dictates of buoyancy scaling. Use of non- 
separable and anisotropic parametric spectral rep- 
resentations can be found in Miiller et al. [1978]. 

Functions A(x = ^^) and B{a) should be nor- 
malized to unity, so that the integrals 

mc/m, oo 

A(x)dx^/^(x)dx = l, 

mi/2m, 



and 



TV 



f B{a)da^ f B{a)da = 1 



f f 

are dimensionless and in practice the limits of inte- 
gration (mi represents vertical mode-1, rric a high 
wavenumber cutoff of approximately 0.1 cpm) are 
approximated for analytic convenience. The high 
wavenumber cutoff, defined by 

2m^Ek{m)dm = —W (17) 

10 
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represents an observed transition to steeper spec- 
tral slopes and likely signifies a dynamical transi- 
tion to a more strongly nonlinear, yet not turbu- 
lent, regime. 

The GM model evolved over time, resulting in 
three major versions, denoted GM72, GM75 and 
GM76: 



GM72 : 



m 



GM75, 76 : 

m 

with 



m 




-1 



if m < m* 
if 771 < m* 



(18) 



M: 



«r(f) 



m^ 



[S, I, J=f 



tzl 
s 



) [1 + (m/m*)^]-*/^ 



(1,2.5,6) 

(2,2,3) 



GM75 
GM76. 



and 



B{a) oc fT-''+2'?(a2 - fy 



(19) 



The gamma function (F) appears in the normaliza- 
tion for the vertical spectrum. All versions of GM 
utilize r = 1 and q = 1/2, so that 



B{a;r = l,q = l/2) 



2/ 



1 



vr ay^ 



P 



Both functions A{x = ^^) and B{a) behave as 
power-laws at high wavenumber and frequency; B 
has an integrable singularity at the inertial fre- 
quency / (which constrains g < 1), and A has a 
plateau for small value of its argument. As ex- 
plained below these features were found to be rep- 
resentative of the oceanic internal wave field. 
The total energy in the model. 



e: 



GM 



30 X 10"^m2s"^ 



(20) 



is based upon fits to observed high frequency spec- 
tra rather than estimates of total observed internal 
wave band energy. 

2.3. The GM model's three incarnations 

The GM model evolved over time as more and 
better data became available. This section docu- 
ments reasons behind the choice of specific param- 
eters. 
2.3.1. Frequency Domain 

While Garrett and Munk [1972] noted variabil- 
ity in the tidal peak, the tides are not part of 
the GM model. The specific shape of the iner- 
tial peak was simply chosen out of analytic conve- 
nience. Variability in the frequency domain power 
law was apparent even then, with a noted alter- 
native being a~^'^ [Webster, 1969]. With many 
possible 'noise' sources such as mooring motion 
and finestructure contamination serving to whiten 
the frequency spectrum, the judicious choice for 



a model of the background spectrum is a redder 
o"~^ (Section 3.3.1). This choice remains stable 
throughout the various model versions. 
2.3.2. Vertical Wavenumber Domain 

The spatial domain information available in the 
late-1960's was limited to stationary current me- 
ters and horizontally towed thermistor chains. The 
direct information available from current meters 
was through spatial lag coherences^. The towed 
data return both a direct estimate of the horizon- 
tal spectrum and vertical lag coherence informa- 
tion. These data supported only a crude repre- 
sentation of the vertical wavenumber spectrum as 
a top-hat model with about 20 equivalent modes 
excited (j* = 20 in GM72). 

The quality of the information improved dra- 
matically in the early 1970's with the introduc- 
tion of near-continuous vertically profiling instru- 
mentation such as the Neil Brown CTD [Brown, 
1974] and Tom Sanford's electric field sensing de- 
vice [Sanford, 1975]. The choice of which data to 
rely upon, though, is quite subjective and the per- 
ception is that finestructure contamination of the 
density field is more problematic than contamina- 
tion of the velocity field, e.g., Polzin et al. [2003]. 
See Section 3.3.2 for a discussion of the Sanford 
[1975] data set. Extant temperature gradient spec- 
tra tended to be white at an intermediate range 
of vertical wavenumbers (roughly 0.01 < m < 0.2 
cpm) and roll off thereafter. GM75, however, relied 
upon even higher wavenumber CTD data from Mil- 
lard [1972], ignored possible inconsistencies with 
CTD data in Hayes [1975], and heavily weighted 
the velocity profile data in Sanford [1975] to argue 
for a high wavenumber slope t = 5/2. Estimates 
of isotherm fluctuations from a vertically profiling 
fioat [Cairns, 1975] suggested a much lower band- 
width (j* = 6) than in the GM72 model. 

GM75 was soon replaced by a more refined 
model. Further analysis of float data by Cairns 
and Williams [1976] suggested revision of the high 
wavenumber power law to m~'^ and a change to 
the functional form l/{m1 + m?). The rationale 
given for choosing such a functional form is sim- 
ply its analytic convenience: l/{rnl + m?) has a 
simple cosine transform, see also Desaubies [1976]. 
Munk [1981] follows suit. This revision is labeled 
GM76. The distinction between l/(m^, + tti)^ and 
\/{ml + m^) is in the rapidity of the roll-off from 
the low-wavenumber plateau region to the high- 
wavenumber asymptote. GM76 has a more rapid 
transition than GM75 for equivalent power laws. 

3. OBSERVATIONS OF INTERNAL WAVE 
SPECTRA 
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3.1. Observational frame'work/preliminary also allow for variable spectral levels with the fo- 
notes 

Observations of the oceanic internal wave field 

are analyzed in this section in order to define the 
extent to which the deep ocean internal wavefield 
is indeed universal. We follow the lead of GM by 
characterizing observed spectra with the paramet- 
ric power law representations (18) and (19) and the spectral parameters. 



cus upon identifying spatial/temporal variability of 
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Figure 1. Geographic depictions of the data sets referred to in this study. Surveys for the Mid-Ocean 
Dynamics Experiment (MODE) and the North Atlantic Tracer Release Experiment (NATRE) are en- 
closed by the circle and square, respectively. This study uses current meter data from near the center 
of these survey regions. The alphabetical keys are: D (Site-D), I (Internal Wave Experiment), L (the 
Long-Term Upper-ocean Study), F (Frontal Air-Sea Interaction Experiment), III-C (PolyMode IIIc) and 
SFTRE (the Salt Finger Tracer Release Experiment) in the North Atlantic; S (SWAPP), P (Patchex), 
M (Misery 1 and 3 as referred to in Cairns and Williams [1976]), C Cairns [1975], E (The offshore ar- 
ray from the Eastern Boundary Currents field program, Chereskin et al. [2000]), MATE (the MidOcean 
Acoustic Transmission Experiment, Levine et al. [1986]) and OS (Ocean Storms, D'Asaro et al. [1995]). 
Bathymetric contours are every 1000 m, alternately bold and thin. The Gulf Stream lies shoreward of 
the 800 m isobath south of Cape Hatteras, cutting between Site-D and the Lotus region. 
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Variability in the frequency domain is quantified 
by fitting 

E{a) = EB{a; r, q = 1/2) oc eoCJ-'+^a^ - f)-^'\ 

(21) 

to the high frequency portion of the observed spec- 
tra. 'High' frequency refers to periods less than 10 
hours, thereby eliminating the inertial peak and 
semi-diurnal tides from consideration. The am- 
plitude Co is reported as a fraction of the energy 
associated with the GM spectrum: 

_ jfEBia;r,q=l/2)da 
^o- EGMN/No ■ ^ ' 

An additional reported statistic is the total in- 
ternal wave band energy in the observed spectrum. 
Internal wave band energy is estimated by integrat- 
ing the observed frequency spectra over frequencies 
greater than approximately 0.7/, thereby including 
broadening of the inertial cusp. 

Variability in the vertical wavenumber domain 
is quantified by fitting variants of (18) to the ob- 
served spectra. The resulting statistics are less sat- 
isfactory than those in the frequency domain as 
low wavenumbers are not well resolved in most in- 
stances. 

3.2. Some w^ords of caution 

Power laws are an attractive and simple way to 
describe complex systems in general and internal 
waves in particular. Some details of the observed 
spectra may defy this simple approach. Our view is 
to treat the parametric representation as a general 
tendency to help summarize and interpret other- 
wise complicated data. With that caveat, there is 
no shortage of reasons for the observed variability: 
nonlinearity, instrumentation, the buoyancy pro- 
file, low frequency flows, tides, bathymetry, hori- 
zontal inhomogeneity, vertical asymmetry and con- 
tamination by either quasi-permanent finestruc- 
ture or self advection (Doppler shifting) within the 
internal wavefield: 

• It matters how you look at the ocean. Instru- 
ment response and processing methods may affect 
the interpretation. An appendix contains remarks 
to aid the reader in negotiating this maze. 

• All that wiggles is not necessarily an in- 
ternal wave. Internal waves typically dominate 
the small scale variability and geostrophically bal- 
anced motions dominate the largest scales. But 
quasi-permanent finestructure has been identified 
as making an increasing contribution to the total 
spectrum at high wavenumber (e.g., Miiller et al. 
[1978]; Polzin et al. [2003], Section 2.1.4). 

• In many cases, the quoted spectral power 
laws are derived from one-dimensional (vertical 
wavenumber and frequency) spectra. Thus there 
is an implicit assumption that the spectra are sep- 



arable, (15). Since both 1-d spectra are red, fre- 
quency spectra are typically dominated by low ver- 
tical wavenumber motions and vertical wavenum- 
ber spectra are dominated by low frequencies. The 
one-dimensional power laws are not identical to the 
2-d power laws if the spectra are non-separable. 

• The sampling strategy matters. Vertical 
wavenumber-frequency domain information are 
seldom obtained from a single instrument, nor are 
the measurements taken instantaneously. Vertical 
profile data obtained as a spatial survey can be 
used to define a spatial mean and residuals, and the 
internal wave contribution identified as the residual 
field. However, this will alias some low-frequency 
variability into the residual fields. Vertical profile 
data obtained as a time series can effectively sep- 
arate the internal wave variability from lower fre- 
quency motions. But the resulting residuals may 
not be representative of the background wavefield: 
that background may be temporarily distorted, for 
example, by wave-mean interactions. 

• Anisotropic propagation in a region of vari- 
able stratification or into/away from a critical layer 
can imply significant transport of energy in the ver- 
tical wavenumber domain and thereby create ap- 
parent structure unrelated to nonlinearity. 

• The results presented here are primarily from 
analyses in depth or pressure coordinates, rather 
than in an isopycnal coordinate system. High fre- 
quency Eulerian spectra are prone to Doppler shift 
contamination at high vertical wavenumber. Use of 
an isopycnal coordinate system limits the Doppler 
shift contamination [Sherman and Pinkel, 1991]. 

Despite these limitations, we believe that the 
patterns of spatial and temporal variability dis- 
cussed below reflect physical reality and are not 
an artifact of processing or interpretation. 

3.3. A Regional Catalog 

Variability of the oceanic internal wave spectrum 
will likely exist on multiple temporal and spatial 
scales, e.g., those scales associated with individual 
wave packets and isolated forcing events. In the 
present study we assume that differences in the 
background spectrum are a product of regional 
differences in forcing and boundary conditions av- 
eraged over the characteristic time scale to dissi- 
pate the energy resident in the internal wavefleld 
(about 50 days) and averaged over the characteris- 
tic spatial scale that internal waves can propagate 
in that time span. The forcing mechanisms and 
propagation scales are, in general, not regionally 
well defined. 

Thus multiple data sets from Site-D (north of 
the Gulf Stream), the Sargasso Sea and the Cali- 
fornia Current System are grouped together (Fig. 
1). In terms of low frequency variability, Site-D ex- 
hibits the largest eddy energy, e.g., Wunsch [1997]. 
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However, variability is apparent in spectral ampli- 
tude. 
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Figure 5. Site-D frequency spectra of horizontal ki- 
netic energy (blue lines). These are the Site-D data that 
appeared in the original GM72 paper. Black curves rep- 
resent fits of (21) with r = 2. The thick vertical lines 
represent the buoyancy frequency cut-off. The spectra 
have been offset by one decade for clarity. 



Eddy energies in the Sargasso Sea and California 
Current System exhibit spatial trends but are com- 
parable to each other. Short narratives are used to 
document regional characterizations of eddy am- 
plitudes, tides and bottom topography. 
3.3.1. Site D - m-^° and o-^-o. 

Site-D was fourth of a series of stations on a sec- 
tion from Woods Hole to Bermuda. Historically, 
the technology for long term current meter deploy- 
ments was developed in Woods Hole, and Site-D 
served as the test bed for much of that instru- 
mentation. Much of the data in the original Gar- 
rett and Munk paper (GM72) came from Site-D, if 
for no other reason than not much else was avail- 
able. Site-D is characterized by a relatively small 
amplitude tide and large low-frequency (predom- 
inantly Topographic Rossby Wave) activity, e.g., 
Hogg [1981]. Site-D is situated under the East- 
ern Seaboard storm track, resulting in an enhanced 
potential for the resonant coupling between wind- 
stress and mixed layer at near-inertial frequencies 
[D'Asaro, 1985]. Site-D may be relatively unique 
because of its geographic location. 

It is apparent that, although many data were un- 
der consideration in Garrett and Munk [1972], the 
GM72 model is a curve fit to the original Site-D 
data presented in Foffonoff [1969] in both ampli- 
tude and 0-^2 power law (Fig. 5)^. The GM72 
model is a good characterization of the spectral 
shape at Site-D throughout the decades (Fig. 2). 
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Figure 2. Frequency spectra of horizontal kinetic en- 
ergy and potential energy (blue and green lines) from 
nearby Site-D. Black curves represent fits of (21) with 
r = 2. Thick vertical lines represent the buoyancy fre- 
quency cut-off. The shape of the fit is stable through- 
out the three decades separating the data in the original 
GM72 work. Note that the amplitude is somewhat lower. 
The Primer-1 spectra have been offset by one decade for 
clarity. 



The original (1967) Site-D data are character- 
ized by larger spectral levels than the other, longer 
term, estimates. Briscoe and Weller [1984] doc- 
ument an annual cycle in high frequency inter- 
nal wave energy using data from the Lotus region 
south of the Gulf Stream, with maximum values in 
winter. This annual cycle is also apparent in the 
Primer data from Site-D (Fig. 3). A seasonal cycle 
in near-inertial energy for this location is discussed 
in detail in Silverthorne and Toole [2009]. 

With the CM spectrum being a curve fit to the 
early (1967) data, and with those data being ob- 
tained in February-March, the GM model is essen- 
tially a description of the winter time spectrum at 
Site-D. The summer time spectrum has lower am- 
plitude with little change in spectral slope at high 
frequency (Fig. 4). We anticipate results presented 
below by stating here that annual average spectral 
levels at Site-D tend to be larger than other places 
and thus the Garrett and Munk model is a poor 
description of the background internal wavefield in 
much of the World Ocean. 
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Figure 3. Site-D time series of horizontal kinetic energy. 
The sohd hne represents 10 times the high frequency en- 
ergy. The dashed line represents the entire internal wave 
band energy estimate. The last 150 days of this record 
were not included in this analysis as they were contami- 
nated by biological fouling of the rotor. 
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Figure 4. Frequency spectra of horizontal kinetic and 
potential energy (blue and green hues) from nearby Site- 
D. The 2-year Primer data set has been divided into high 
and low amphtude states. The low amplitude state has 
been offset one decade for clarity. Black curves repre- 
sent fits of (21) with r = 1.90 (high amplitude state) and 
r = 1.95 (low amplitude state). The thick vertical line 
represents the buoyancy frequency cut-off. 



Despite the vertical spectrum being defined us- 
ing data obtained elsewhere, recent vertical profile 
data data from Site-D are remarkably consistent 
with the GM76 model [l/(m2-hm2)] (Fig. 6). Win- 
tertime conditions exhibit both enhanced spectral 
levels (amplitude factors of 2.75 vs 1.75) and rel- 
atively more variance at low modes (j* =4 — 5 
vs j* = 10) than summertime data. Winter time 
conditions also exhibit larger ratios of kinetic-to- 



potential energy at high wavenumber, implying an 
increased input of near-inertial energy during win- 
ter time and relaxation to higher frequencies. 
3.3.2. The Sargasso Sea 

A large number of experiments have been lo- 
cated in the Sargasso Sea over the Hatteras Abyssal 
Plain. On the southern side of the Gulf Stream, 
this region exhibits an energetic eddy field having 
significant north-south gradients. Eddy energy lev- 
els are typically less than noted at Site-D. A tidal 
(M2) peak is apparent in the temperature and ve- 
locity spectra. Miiller et al. [1978] find that fluctu- 
ations at this frequency have larger characteristic 
vertical scales than the internal wave continuum, 
and there is evidence of similar features at the first 
several harmonics. From current meter data at 
28°N, 70°W Noble [1975] and Hendry [1977] es- 
timate net fluxes at M2 to be to the southeast and 
infer the source to be the Blake Escarpment, near 
the western boundary. Alford and Zhao [2007], on 
the other hand, document net semidiurnal fluxes 
to the north-northwest (at 31°N, 69°30'W ) and 
southwest (at 34°N, 70° W ). 

The bottom near mid-basin is well sedimented 
and smooth at 28°N, 70°W, the locus of the 
Mid-Ocean Dynamics Experiment and the Internal 
Wave Experiment. Rougher topography is noted 
to the east^. The Poly Mode Local Dynamics Ex- 
periment current meter array (at 31°N, 69°30'W 
) is likely situated above relatively rough terrain. 
Bathymetry is relatively smooth at 34°N, 70°W, 
the locus of the Long Term Upper Ocean Study. 

The buoyancy frequency profile has a relative 
minimum in N{z) associated with the 18° water 
thermostad at about 300 m water depth. The main 
thermocline exhibits nearly constant stratification 
between 500-1000 m. 

MODE 

Conducted during March-July of 1973, the Mid- 
Ocean Dynamics Experiment (MODE) was one of 
the first concentrated studies of mesoscale ocean 
variability. The experiment featured arrays of 
moored current meters, neutrally buoyant floats, 
standard hydrographic station techniques and the 
use of novel vertically profiling instrumentation. 
An extensive array of current meter moorings was 
deployed in a 300 km radius centered about (28° 
N, 69°40' W). Data return from the current meters 
was Umited [Mode Atlas, 1977]. 

While designed primarily to investigate low fre- 
quency motions, the experiment returned a great 
deal of information about internal waves. Verti- 
cal profiles of horizontal velocity obtained during 
May and June with a free-falling instrument us- 
ing a electric field sensing technique Sanford [1975] 
provided, for the first time, direct estimates of the 
high vertical wavenumber structure of the ocean 
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Figure 6. Site-D vertical wavenumber spectra of horizontal kinetic and potential energy. The black line 
represents a fit of the GM76 spectrum. Spectral estimates at the lowest wavenumber (enclosed circles) 
were made using a modal fit. These data were recently acquired as part of a long term climate monitor- 
ing project, Station-W (unpublished data from J. Toole, personal communication 2004). Velocity and 
density profiles were obtained with the Moored Profiler. Information regarding the internal wavefield 
is returned by burst sampling 4 times using a 9.5 hour sampling interval, then waiting 5 days before 
repeating. The spectrum presented here represents departures from the burst means. See the appendix 
for further details. Departure from the curve fit at vertical wavelengths of 10 m and smaller is interpreted 
as noise. 



internal wavefield. These data are dominated by 
near-inertial frequencies, and a rotary decomposi- 
tion in the vertical wavenumber domain Leaman 
and Sanford [1975] reveals a large excess of clock- 
wise {cw) phase rotation with depth. Clockwise 
phase rotation with depth is a signature of down- 
ward energy propagation for near-inertial waves. 

Despite the evidence of excess downward energy 
propagation that was interpreted in terms of atmo- 
spheric generation [Leaman, 1976], these data were 
assumed to be representative of the background in- 
ternal wavefield. A subset obtained as a time series 
of 4^ days provided the basis for a revision to the 
isotropic Garrett and Munk spectral model, GM75, 
with high wavenumber asymptote of m~^''^ . 

The MODE profile data were included in our 
original study [Lvov et al., 2004] reporting the co- 
variability of frequency- wavenumber power laws. 
Further investigation suggests the characterization 
of these data as representing the background wave- 
field is problematic. Several coherent wavepackets 
dominate the high wavenumber energy content and 
these packets can be interpreted as being strained 
by the mesoscale eddy field and propagating into 
a critical layer [Polzin, 2008a]. Thus wave-mean 
interactions acting in conjunction with the down- 
ward propagating waves may be responsible for 
creating structure in the vertical wavenumber do- 



main uncharacteristic of the background spectrum. 

IWEX - /fc-2.4±0.4 ^^^ ^-1.75 

The Internal Wave Experiment (IWEX) rep- 
resents an early attempt to estimate a vector 
wavenumber-frequency spectrum with a minimum 
of assumptions. Current and temperature data 
were obtained with 17 Vector Averaging Current 
Meters (VACMs) and 3 Geodyne 850 current me- 
ters from a taut, three-dimensional, trimoored ar- 
ray. 

Miiller et al. [1978] assumed the spectrum could 
be factorized as: 

E'^ik, if, a) = E''{a)A{k; a)S{ip; a) 

with (T = lb denoting the sign of vertical en- 
ergy propagation, A the normalized horizontal 
wavenumber (k) distribution and S a normal- 
ized azimuthal {(p) distribution. Their horizontal 
wavenumber distribution was assumed to have a 
parametric representation of: 



A{k;a)oc [l + ( 



rC rC'Fi 



T\s]-t/s 



f^^ 



for k > kp and ^ = otherwise. The parameters 
are horizontal wavenumber scale A;*, low wavenum- 
ber cut-off kp, shape factor s, and high wavenum- 
ber slope t. An inverse analysis of the spatial lag 
cross spectra (the lag-coherence is proportional to 
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Figure 9. Iwex frequency spectra of horizontal kinetic 
energy and potential energy (blue and green lines). Black 
curves represent fits of (21) with r = 1.75. The thick ver- 
tical line represents the buoyancy frequency cut-off. 



mode. At near-inertial frequencies kp corresponds 
to the third vertical mode, indicating significantly 
reduced near-inertial energy in the first two modes. 
This pattern is consistent with the vertical profile 
data presented here which depict ratios of Ej. to Ep 
at 100 (1000) m. vertical wavelengths that are sig- 
nificantly larger (smaller) than in the GM model, 
for which E^/Ep = 3. 

Miiller et al. [1978] avoid interpreting the IWEX 
spectrum in the context of how it might relate to 
the background oceanic spectrum. This is not sur- 
prising given the short duration of the data set and 
the relative lack of quality data sets at that time. 
Here, however, we promote the IWEX spectrum as 
a more realistic representation of the Sargasso Sea 
background spectrum than the GM model. 



the Fourier transform of the energy spectrum) was 
performed to estimate the various parameters in 
their proposed spectrum. In particular, Miiller 
et al. [1978] find that the horizontal wavenum- 
ber energy spectrum depends upon wavenumber 
as /c~^'^ , independent of frequency. For contin- 
uum (/^ <^ (J^ <^ N'^) frequencies the horizontal 
and vertical wavenumber power laws are identical. 
Thus E{m) oc m"^'^^'''^. The observed frequency 
spectra of both velocity and temperature are char- 
acterized by the power law a^^'"^^, Briscoe [1975b] 
(Fig. 9). 

The Miiller et al. [1978] analysis includes an es- 
timate of finestructure contamination. That is, 
their horizontal wavenumber spectrum includes 
only contributions which are consistent with lin- 
ear internal waves. An increasing contamination 
with increasing frequency (decreasing horizontal 
scale) is apparent. As they make no distinc- 
tion between permanent finestructure contamina- 
tion and Doppler shifted internal waves, much of 
their finestructure contamination could be internal 
waves [e.g., Sherman and Pinkel [1991]; Polzin et 
al. [2003]]. Thus the horizontal wavenumber do- 
main power law estimate may be biased too steep. 

Several features of the IWEX data set reoccur 
in many of the other data sets analyzed here. The 
first is that estimates of amplitude based upon fits 
to the high frequency part of the spectrum tend to 
be significantly lower than the GM model whereas 
estimates of total energy tend to be more nearly 
equal to the total energy in the GM model. The 
inference is that the inertial cusp of the GM model 
is not sufficiently strong to describe the oceanic 
peak. While this has been noted before (e.g., Fu 
[1981]), the IWEX analysis indicates the excess 
near-inertial energy has smaller characteristic ver- 
tical scales: At super-tidal frequencies the peak 
wavenumber kp corresponds to the first vertical 
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LOTUS 

One goal of the Long-Term Upper-Ocean Study 
(LOTUS) was to document variability in the in- 
ternal wavefield over several seasons in order to 
investigate its association with a variety of forcing 
mechanisms and environmental conditions. Mete- 
orological data were obtained with a surface moor- 
ing and buoy. Only data from subsurface moorings 
are described here. The region is notable for ex- 
treme air-sea buoyancy exchange resulting in the 
production of Eighteen Degree Water, e.g., Kwon 
and Riser [2004] 

Briscoe and Weller [1984] document an annual 
cycle of high frequency internal wave energy at 200- 
500 m depths from the first year's data. The en- 
ergy was noted to vary from a half or a third to 2 
or three times the mean at each depth. That sig- 
nal is repeated in the second year's data (Fig 7). 
Apart from the seasonal cycle, internal wave en- 
ergy varies on time scales of several weeks. The re- 
lation between this short time scale variability and 
low frequency (sub-inertial) shear has been noted 
since MODE (see Ruddick and Joyce [1979] and 
references therein). 

In this study of how nonlinearity may shape and 
form the internal wave spectrum, there is a po- 
tential link between spectral amplitude and spec- 
tral shape. The Lotus spectra in the 325-350 m 
depth range were averaged over the time periods 
indicated in Fig. 7. Spectra from the more en- 
ergetic time periods are slightly less steep {a~^''^^ 
vs. a~^'^^) for a factor of 2.5 difference in spectral 
level. If there is a dynamical link between spectral 
level and spectral shape, that link is subtle relative 
to the gamut of variability discussed here. 
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Figure 7. Lotus time series of high-frequency energy 
from 325 (year-1) and 350 (year-2) m water depth. This 
depth range is occupied by 18 Degree Water and repre- 
sents a local minimum in the stratification rate. Mixed 
layer depths were observed to be smaller than 300 m over 
the duration of the data set. The solid Une represents 10 
times the high frequency energy. The dashed line rep- 
resents the entire internal wave band energy estimate. 
The horizontal lines indicate the time periods over which 
spectra were averaged into bins of high and low energy 
states of Fig. 8. 
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Figure 8. Lotus frequency spectra of horizontal kinetic 
energy and potential energy (blue and green lines) for 
the high and low energy states depicted in Fig. 7. Black 
curves represent fits of (21) with r = 1.75 and r — 1.85. 
The thick vertical line represents the buoyancy frequency 
cut-off. Temporal variability is dominated by variability 
in the overall amplitude of the spectra rather than the 
shape (power law). 



FASINEX - m-2-3 and o-^-^s 

The Frontal Air-Sea Interaction Experiment 
(FASINEX) was designed to investigate the re- 
sponse of the upper ocean to atmospheric forc- 
ing in the presence of oceanic fronts. An array 
of surface and subsurface moorings with VMCMs, 
VACMs and Profihng Current Meters (PCMs) was 
deployed in the subtropical convergence zone of 
the Northwest Atlantic ( approximately 27° N, 70° 



W) from January to June of 1986, Weher [1991] 
and Eriksen et al. [1991]. The moored data in 
this study are taken from a long-term, subsurface 
mooring at 28° N. These data document an annual 
cycle in internal wave energy that is substantially 
reduced from the Lotus time series 500 km north 
(Figure 10). Frequency spectra (Fig. 11) at depths 
of 556 and 631 m are defined by a power law a~'^ of 
r = 1.85 and amplitudes significantly smaller than 
the GM model. 

Vertical profiles of horizontal velocity and den- 
sity were obtained during February-March using 
the High Resolution Profiler (HRP), Polzin et al. 
[1996]. The vertical profiles, obtained primarily 1° 
north of the moored array as part of a spatial sur- 
vey, revealed a complex pattern of variability as- 
sociated with the frontal velocity structure in the 
upper 250 m. Here we report results concerning 
data from depths of 250-1000 m. Those data are 
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Figure 10. FASINEX time series of high-frequency en- 
ergy from 556 and 631 m water depth. The solid line rep- 
resents 10 times the high frequency energy. The dashed 
line represents the entire internal wave band energy esti- 
mate. 
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Figure 11. FASINEX frequency spectra (blue and green 
lines) of horizontal kinetic energy and potential energy 
from the main thermocline (600 m). Black curves rep- 
resent fits of (21) with r — 1.75. The thick vertical line 
represents the buoyancy frequency cut-off. 
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fit with a spectrum having an asymptotic roll-off 
of m-2-3 (Fig. 13). 

Weller et al. [1991a] document a slight excess of 
clockwise phase rotation with depth over counter- 
clockwise shear variance, a signature of excess 
downward propagating near-inertial energy. The 
quoted power law in the vertical wavenumber do- 
main is at a depth for which N"^ varies by less than 
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Figure 13. FASINEX-vertical wavenumber spectra of 
horizontal Icinetic and potential energy, A^-scaled and 
stretched under the WKB approximation to No = 3 
cph. These data were obtained 100-200 km north of the 
FASINEX moored array. Sampling was intermittent in 
both space and time. Note that the low- wavenumber 
spectral estimates have typically smaller ratios of -B^ and 
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Figure 14. PolyMode IIIc frequency spectra (blue and 
green lines) of horizontal kinetic energy and potential en- 
ergy from the main thermocline (300 m). Black curves 
represent fits of (21) with r — 1.90. The thick vertical 
line represents the buoyancy frequency cut-off. 



a factor of 2, so that, biases associated with linear 
wave propagation in nonuniform A^ are likely small. 
The instrument does not, however, return robust 
velocity estimates at large vertical scales and hence 
diagnosing interactions with the thermocline scale 
geostrophic shear is problematic. 
3.3.3. The wrestern Tropical Atlantic 
SFTRE and PolyMode IIIc - m'^ "^ and a"^^ 

Extant frequency domain data from thermocline 
regions in the western Tropical Atlantic are limited 
to those obtained as part of the PolyMode pro- 
gram. The motivation for PolyMode Array III was 
exploration of low frequency variability in what 
were perceived as dynamically distinct regions of 
the North Atlantic. Cluster IIIc was placed in the 
eastward flowing North Equatorial Current with 
the intent of examining the low frequency variabil- 
ity for characteristics of eddy generation by baro- 
clinic instability, Keffer [1983]. The cluster is sit- 
uated over the northwestern extension of the De- 
merara Abyssal Plain. Locally the bottom is well 
sedimented and relatively flat, though rough to- 
pography lies immediately to the east. Data pre- 
sented here are an average of the three southeast- 
ern moorings (80-81-82). 

Frequency spectra (Fig. 14) roll-off less steeply 
than (T~^. The simple power law characterization 
(19) overestimates the observed spectral density at 
frequencies smaller than semi-diurnal. An annual 
cycle in either high frequency or internal wave band 
energy is not apparent in the time series (Fig. 15). 

Vertical proflle (HRP) data used here were ob- 
tained during Nov. 2001 as part of a sampling 
survey for the Salt Finger Tracer Release Exper- 
iment (SFTRE). Warm, high-salinity Subtropical 
Underwater overlies cooler, fresher Antarctic In- 
termediate Water [Schmitt et al., 2005]. The situa- 
tion is unstable to the salt-fingering form of double 
diffusive instability and a staircase layering of the 
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Figure 15. PolyMode IIIc time series of high-frequency 
and internal wave-band energy from nominal 300 m wa- 
ter depth. The solid line represents 10 times the high 
frequency energy. The dashed line represents the entire 
internal wave band energy estimate. 
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temperature, salinity and density profiles is present 
over much of the region. The profiles examined 
here were collected between (14 — 16° N, 50 — 57° 
W) on the northern edge of the survey region. 

Vertical wavenumber E^ spectra (Fig. 12) roll- 
off more steeply than m~^ at high wavenumber. 
The spectra are anomalous in that potential energy 
estimates exceed kinetic energy at high wavenum- 
ber, a feature that we attribute to the staircase 
features. The parametric fit (18) produces a low- 
wavenumber roll-off equivalent to mode-9. Esti- 
mates of Ep dominate E/^ at low wavenumber. 
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Figure 12. SFtre2-vertical wavenumber spectra of 
horizontal kinetic and potential energy, A^-scaled and 
stretched under the WKB approximation to No — 3 
cph. These data were obtained on the northern part 
of the year-1 sampling survey grid. Note that the low- 
wavenumber spectral estimates have typically smaller ra- 
tios of Ek and Ep. 



3.3.4. The southeast Subtropical North At- 
lantic 

NATRE and SUBDUCTION - m^^^ or 
m"^-'^^ and a~^-^^ 

Vertical profile (HRP) data were collected dur- 
ing April 1992 southwest of the Canary Islands as 
part of an initial site survey for the North Atlantic 
Tracer Release Experiment (NATRE) . The bulk of 
the data were obtained as part of a 400 x 400 km 
grid centered about (26°N, 29° W). Apart from a 
minor seamount with a summit at 3000 m water 
depth, the bottom is relatively featureless near the 
survey domain. The region is further notable for 



the production of a relatively minor water mass, 
Madeira Mode Water [Weller et al., 2004]. 

A rotary decomposition of the vertical profile 
data does not return a consistent pattern of phase 
rotation with depth as a function of wavenumber. 
Thus, despite the fact that the buoyancy profile in 
the southeast part of the subtropical gyre decreases 
monotonically from the mixed layer base through 
the main thermocline, a possible bias of the spec- 
tral slope associated with a purely linear response 
to wave propagation in variable N{z) is unlikely. 

The vertical wavenumber kinetic energy spec- 
tra (Fig. 16) are white at low wavenumber and 
roll-off more steeply than m~^ for vertical wave- 
lengths smaller than 200 m {E^ oc m"^'^^). The 
observed potential energy spectra exhibit three 
salient features in contrast with the kinetic energy 
spectrum. First, kinetic and potential energy are 
nearly equal at the largest resolved vertical wave- 
lengths, about 1000 meters. The spatial survey is 
of sufficient lateral extent that, even though low 
frequency variability is dominated by barotropic 
and low mode (mode-1 and mode-2) contributions 
(Wunsch [1997]), low frequency variability may 
contribute to the observed spectra at the largest 
resolved scales (equivalent to about mode-3 in Fig. 
16). Analysis of a subset of these data obtained on 
a grid with 1 nautical mile spacing returns a con- 
sistent signature: Ek/Ep ratios of 2-3 are found at 
the lowest two resolved vertical wavelengths. Sec- 
ond, the spectra diverge so that Ek is about an 
order of magnitude larger than Ep at 100 m ver- 
tical wavelengths. Third, Ep does not roll-off as 
quickly at high wavenumber, as can be inferred 
from increasing ratios between Ep and Ej. at ver- 
tical wavelengths smaller than 10 m. Polzin et al. 
[2003] argue that these increasing ratios at small 
scales are largely associated with an increasing con- 
tribution of quasi-permanent finestructure (Section 
2.1.4). When the quasi-permanent contribution is 
subtracted from the observations, both Ek and Ep 
roll-off at about the same rate, E^ (x Ep cc m~'^''^°. 
The low wavenumber excess of Ep remains. The 
vertical wavenumber spectrum indicates a non- 
separable spectrum, with excess near-inertial con- 
tent at high wavenumber. 

The vertical profile data were obtained from the 
vicinity of the center mooring of the Subduction ar- 
ray, Weller et al. [2004]. VMCMs deployed on sur- 
face moorings document obvious inertial and tidal 
peaks, and peaks at several harmonics (Fig. 19). 
The frequency spectrum is not succinctly charac- 
terized in terms of a single power law. At low fre- 
quencies (1 < o" < 6 cpd) the spectrum exhibits 
a plateau (ignoring the tidal peak and harmonics) 
and at higher frequencies (6 < a < 48 cpd) rolls 
off as cr~^'^^. Energy ratios [Ek{(7) / Ep{a)] agree 
with linear kinematics (9), Internal wave band ra- 
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Figure 16. Natre vertical wavenumber kinetic Ek and potential Ep energy spectra. These 100 profiles 
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tios [J Ek{(7)da/ J Ep{a)da] are between 2-3, con- 
sistent with a red vertical wavenumber spectrum 
having httle near-inertial content at low wavenum- 
ber. Unlike the Sargasso Sea, no seasonal cycle is 
apparent in either internal wave energy time series 
(Fig. 18). 

NATRE is special among these data sets in that 
it exhibits a strikingly large M2 tide and the re- 
gion is a relative minimum in eddy energy, Wun- 
sch [1997]. An obvious hypothesis is that structure 
in the spectral domains reflects the detailed path- 
way energy takes in draining out of the tidal peak. 
A second hypothesis is suggested by noting that 
the transition between the two regimes (roughly 
6 cpd) corresponds to the stratification rate near 
the bottom boundary (A^ = 3 cpd). Waves of 



higher frequency never reach the bottom and will 
not be attenuated by viscous processes in the bot- 
tom boundary layer. 

3.3.5. The southeast Subtropical North Pa- 
cific 

A significant number of data sets have been ob- 
tained offshore southern California. The major 
geographic feature here is the California Current 
System consisting of a deep poleward flow along 
the continental slope, equatorward flow at the shelf 
break and a 600-km wide band of mesoscale vari- 
ability extending from the shelf and fading into 
the northeast Pacific (e.g., Chereskin et al. [2000]). 
There are, as well, gradients in mesoscale variabil- 
ity along the coast, with the most northern of the 
data sets discussed here being closest to a regional 
maximum in surface kinetic energy estimated from 
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surface drifter data Chereskin et al. [2000]. The 
data discussed here are more than 400 km distance 
from the shelf break. The buoyancy frequency 
[A^(z)] typically attains a maximum at 100 m wa- 
ter depth and decreases by about a factor of three 
over 300 m. This represents a large gradient. 

Data from this region include the oscillating 
float estimates of isopycnal displacement provid- 
ing frequency spectra and vertical lag coherence 
estimates reported in Cairns [1975] and Cairns and 
Williams [1976]. These data prompted the revision 
of GM75 model to the GM76 form. See Section 2.3 
for further discussion. 

Here we utilize moored current meter data from 
the Eastern Boundary Currents program and ver- 
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Figure 18. Natre time series of high-frequency and 
internal wave-band energy from nominal 300 m water 
depth. The solid line represents 10 times the high fre- 
quency energy. The dashed hue represents the entire in- 
ternal wave band energy estimate. 
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Figure 19. Natre frequency spectra (thin lines) of hori- 
zontal kinetic energy and potential energy from the main 
thermocline (300 m). Observed spectra (soHd lines) with 
vortical mode spectra (dashed) superimposed. The non- 
propagating vortical fluctuations are assumed to be pas- 
sively advected by the mesoscale field in this representa- 
tion. See Polzin et al. [2003]. Inertial, semi-diurnal and 
harmonic peaks are noted. Black curves represent fits of 
(21) with r = 1.50. The thick vertical hue represents the 
buoyancy frequency cut-off. 



tical profile data from the Patches Experiment 
(PATCHEX; Sherman and Pinkel [1991]) and the 
Surface Wave Process Program (SWAPP; Weller et 
al. [1991b]). PATCHEX was a multi-investigator 
study of the space and time structure of mixing 
events. SWAPP was nominally focussed upon up- 
per ocean processes, but data was obtained from 
thermocline regions as well. These latter two 
field programs featured yo-yo CTD and Doppler 
sonar data sets obtained from the research plat- 
form FLIP. These simultaneous vertical profiles of 
velocity and density provide a unique representa- 
tion of the upper ocean wavefield and permit 2-D 
vertical wavenumber-frequency spectra to be eval- 
uated. This has been done in both depth and den- 
sity coordinates, with the latter being consistent 
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Figure 20. Eastern Boundary Currents time series of 
high-frequency and internal wave-band energy from nom- 
inal 600 m water depth. The solid line represents 10 
times the high frequency energy. The dashed line rep- 
resents the entire internal wave band energy estimate. 
Peaks in high frequency and wave-band energy during 
April and November of 1993 coincide with the presence 
of anti-cyclonic eddies, Chereskin et al. [2000]. 
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Figure 21. Eastern Boundary Currents frequency spec- 
tra (blue and green lines) of horizontal kinetic energy 
and potential energy from the main thermocline (600 
m). Black curves represent fits of (21) with r = 2.00. 
The thick vertical line represents the buoyancy frequency 
cut-off. 
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with and directly comparable to the theoretical 
analysis presented in the companion manuscript 
[Lvov et al., 2010]. The use of this coordinate sys- 
tem will also likely limit contamination by quasi- 
permanent density finestructure. Quantifying this 
is difficult, as is quantifying the effects of Doppler 
shifting by sub-inertial currents. Results from the 
density data are used exclusively here as the ve- 
locity data are subject to contamination from in- 
strumental noise and beam separation effects at 
high wavenumber and frequency, Appendix B2. A 
further cautionary note is that sampling during 
PATCHEX may not have been of sufficient tempo- 
ral or spatial extent for the observations to reflect 
the background wavefield. 

Eastern Boundary Currents - cr ^ 

A coherent moored current meter array was de- 
ployed for a two year period (1992-1994) at the off- 
shore edge of the California Current as part of the 
Eastern Boundary Currents field program [Chere- 
skin et al., 2000]. These data represent the only 
available long-term current meter records in the 
vicinity of the other data from the northeast Pa- 
cific discussed herein. As those data are from short 
term (several week) field programs, it is a natural 
question to ask whether those data are representa- 
tive of the long-term average (Fig. 21). The data 
(both short and long term) share a common fre- 
quency dependence of o""^ at high frequency. Vari- 
ability of both high-passed and wave-band energy 
levels is associated with eddy variability rather 
than a seasonal cycle (Fig. 20). 



val 150-406 m. No information regarding vertical 
symmetry or mean currents is available. 



m 
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PATCHEX^ 

Data from a freely-falling vertical profiler MSP 
were obtained during PATCHEX. Gregg and San- 
ford [1988] note that the velocity profile data are 
remarkable only for the absence of coherent near- 
inertial wave packets. Gregg et al. [1993] document 
a shear spectrum that varies from white at pres- 
sures of 250-550 db to shghtly blue (575-925 db) 
over vertical wavelengths of 10-100 m (Fig. 17). 

PATCHEX^ - m-i-^5 and a'^-^^ - cr-^-O. During 
PATCHEX two CTDs were deployed from the R/P 
FLIP to obtain density profiles over 0-560 m wa- 
ter depth, with a cycle time of 180 s. The density 
data considered here come from the depth inter- 




Figure 17. Patchex vertical wavenumber shear spectra 
using the free-fall profiler MSP. The dashed line repre- 
sents rn? /{mi -\- rn?), with j, = 3 and an m~^ depen- 
dence at vertical wavelengths smaller than 0.1 cpm. The 
domain of interest here is 0.01 < m < 0.1 cpm. (Figure 
extracted from Gregg et al. [1993], their Fig. 6, permis- 
sion from J. Phys. Oceanogr., Copyright 1993 American 
Meteorlogical Society) . 



For vertical wavenumber 1/256 < m < 1/10 
cpm, the density gradient spectrum rolls off as 
m' 



-1.75 



over 0.125 < a < 1.4 cph (Fig. 22) This 
range of wavenumbers and frequencies is charac- 
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■a 



terized by a frequency dependence of a 
The cumulative spectrum (integrated over all fre- 
quencies f < (T < N) is slightly blue for 0.02 < 
rn, < 0.1 cpm. These power law relations were de- 
termined graphically with use of a straight edge. 

SWAPP - m-i-9 and cr-^-O 

During SWAPP a pair of CTDs deployed from 
FLIP returned density profiles from 5 to 420 m 
with a vertical resolution of about 1.5 m and cy- 
cle time of 130 s for a duration of 12 days dur- 
ing March, 1990. The data cited here come from 
depths of 50-306 m. 

For vertical wavenumbers 1/128 < m < 1/10 
cpm, the density gradient spectrum rolls off as a~'^ 
over 0.2 < cr < 4 cph (Fig. 23). For frequen- 
cies of 0.02 < (T < 2 cph, the vertical wavenumber 
spectrum is nearly white, m^-^, so that the energy 
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Figure 22. Patchex 2-D vertical wavenumber - fre- 
quency spectra of isopycnal displacement gradient. The 
density data from which displacement is estimated were 
obtained with a rapid-profiling CTD from the research 
platform FLIP. The quoted power law dependencies have 
been determined graphically. (Figure extracted from 
Sherman and Pinkel [1991], their Fig. 6, permission from 
J. Phys. Oceanogr., Copyright 1991 American Meteoro- 
logical Society.) 

2 *cl 1 (d) 



&io' 







E^nni'lnon-niir 


■ 


|i.'l>il Lii:c;jv 


- 


^ 


j 


\ 


i 


"^'Z^I'' \^ 


. 


iii-iu \n 




■14.3: \ 

31-lf. ^ 


: 


1£-IU 


III 



t3i||ii;>lrh^Tk|i':liii 



10 ItJ 10" 




10 10 



Figure 23. Swapp frequency spectra of isopycnal dis- 
placement gradients. The density data from which dis- 
placements are estimated were obtained with a rapid- 
profiling CTD from the research platform FLIP. The 
quoted power law dependencies have been determined 
graphically. (Figure extracted from Anderson [1992], 
their Fig. 2.9, permission from S. Anderson. Copyright 
1992.) 



spectrum rolls off as m~^'^ (Fig- 24). As with the 
Patchex data from FLIP, these power law relations 
have simply been determined using a x-by-eye pro- 
cedure. 

The buoyancy frequency [A^(2)] attains a maxi- 
mum at 100 m water depth and decreases by about 
a factor of three over 300 m. At periods greater 
than 8 hours, the shear spectra indicate a domi- 
nance of downward energy propagation, Anderson 
[1992]. Thus the vertical spectrum may be influ- 
enced by purely linear kinematics, rather than non- 
linearity. The shear spectrum is not easily charac- 
terized. It exhibits a peak at about 0.03 cpm as- 
sociated with downward propagating near-inertial 
waves, Anderson [1992]. 

3.3.6. The northeast Subpolar North Pa- 
cific 



Figure 24. Swapp vertical wavenumber spectra of isopy- 
cnal displacement gradients. The density data from 
which displacement is estimated were obtained with a 
rapid-profiling CTD from the research platform FLIP. 
The quoted power law dependencies have been deter- 
mined graphically. (Figure extracted from Anderson 
[1992], their Fig. 2.7, permission from S. Anderson. 
Copyright 1992.) 
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MATE - m-^-1 and a' 

An extensive set of moored and profiling (both 
vertical and horizontal) observations was taken 
during June- July of 1977 as part of the Mido- 
cean Acoustic Transmission Experiment (MATE), 
Levine et al. [1986]. The experiment was located 
at 46°46'N, 130°47'W, midway between Cobb 
seamount, which rises to within 30 m of the sur- 
face, and the smaller Corn seamount which rises 
to depths of 1000 m. The seamounts are sepa- 
rated by approximately 18 km and the local wa- 
ter depth is 2200 m. The seamounts may serve 
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Figure 25. MATE frequency spectra of vertical dis- 
placement. Solid and dashed curves are the MATE and 
GM models, respectively. (Figure extracted from Levine 
et al. [1986], their Fig. 5. Permission from J. Geophys. 
Res., , copyright 1986 American Geophysical Union. ) 
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as significant generators of an internal tide. Cobb 
seamount lias been noted to support a trapped di- 
urnal tide [Codiga and Eriksen, 1997; Lueck and 
Mudge, 1997]. The MATE observations focused 
upon subthermocline depths of 900-1300 m away 
from the seamounts. 

Current meter data document frequency domain 
power laws of a~^''^ (Fig.s 25 and 26). Tempera- 
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Figure 26. MATE frequency spectra of horizontal ki- 
netic energy. Solid and dashed curves are the MATE and 
GM models, respectively. (Figure extracted from Levine 
et al. [1986], their Fig. 8. Permission from J. Geophys. 
Res., , copyright 1986 American Geophysical Union. ) 
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Figure 27. MATE vertical wavenumber spectra of ver- 
tical displacement. (Figure extracted from Levine et al. 
[1986], their Fig. 4. Permission from J. Geophys. Res., , 
copyright 1986 American Geophysical Union. ) 



ture profiles from a Bissett-Berman 9040-5 CTDSV 
were used along with a locally tight 9 — S relation to 
estimate vertical displacement. Vertical wavenum- 
ber domain power laws of m"^'^ are quoted for 
6 X 10^3 < m < 6 X 10^2 cpm (Fig. 27). Hor- 
izontal wavenumber spectra of vertical displace- 
ment (again estimated from temperature) obtained 
from the Self-Propelled Underwater Research Ve- 
hicle (SPURV) roll off as | k |~^ for horizontal 
wavenumbers 4.7 x 10"*^ <| A; |< 9.0 x 10~'^ cpm 
(Fig. 28). A combined analysis of spatial lag co- 
herences returns a roll off of j* = 6. 
Strex and Ocean Storms - m ^^ and o"^^^ 

The northeast Subpolar North Pacific has been 
the subject of a number of experiments focussing 
on near-inertial internal wave generation and 
mixed layer processes. The experiments [vertical 
profile data from the Storm Transfer and Response 
Experiment (STREX, D'Asaro [1984]) taken dur- 
ing November of 1980 at 141°W, 50° N and cur- 
rent meter data from the Ocean Storms Experi- 
ment D'Asaro et al. [1995] at 139.25°W, 47.5° N 
are examined here] are located just south of the 
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Figure 28. MATE horizontal wavenumber spectra of 
vertical displacement estimated from SPRUV. Solid and 
dashed curves represent are north-south and west-east 
runs, respectively. (Figure extracted from Levine et al. 
[1986], their Fig. 9. Permission from J. Geophys. Res., , 
copyright 1986 American Geophysical Union. ) 
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climatological storm track in order to emphasize 
the couphng between inertial motions in the mixed 
layer and an inertial time scale anticyclonic rota- 
tion of the windstress vector associated with small 
scale atmospheric disturbances. The region is also 
one of relatively low eddy energy and so minimizes 
the confounding effect of mesoscale eddy - internal 
wave interactions. 

Current meter data from 200 m water depth 
from the Ocean Storms Experiment clearly doc- 
uments the input of near-inertial internal wave 
energy associated with the passage of individual 
storms (Fig. 29). The presence of a seasonal cy- 
cle in inertial energy, though, is obvious only at 
greater depths (500 m). In a broad survey of his- 
torical current meters, Alford and Whitmont [2007] 
find that the seasonal cycle in near-inertial energy 
is most evident at latitudes of 25 — 45°, with de- 
creasing seasonal modulations north of 45° . A pos- 
sible interpretation is a regional emphasis on prop- 
agation. 
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Figure 29. Ocean Storms time series of high-frequency 
and internal wave-band energy from nominal depths of 
200 (upper panel) and 500 m (lower panel) water depth. 
The solid line represents 10 times the high frequency en- 
ergy. The dashed line represents the entire internal wave 
band energy estimate. The passage of storms resulting in 
strong forcing of mixed layer inertial currents is denoted 
by asterisks. 



High frequency energy at 500 m also exhibits a 
seasonal cycle. Unlike seasonal cycles in the west- 
ern North Atlantic, here the high frequency energy 
appears to lag near-inertial energy by several lO's 
of days. A possible interpretation is that the pres- 
ence of such a phase lag is related to the absence 
of an energetic eddy field. 

The frequency spectrum fit indicates a spectrum 
redder than GM. This is paired with a vertical 
wavenumber spectrum of horizontal velocity also 
steeper than GM (Fig. 31). This represents the 
sole example of power laws that do not covary. 
3.3.7. The Arctic 
AIWEX - m-2-25 and a'^-^ 

The Arctic Internal Wave Experiment (AIWEX) 
was conceived as an attempt to study an anoma- 
lous internal wavefield. The experiment took place 
about 350 km north of Prudhoe Bay, Alaska, 
with sampling extending from March to May of 
1985. Total energy levels were a factor of 0.02 to 
0.07 times smaller than in the GM model spec- 
trum. Possible reasons for the low energy levels 
are weaker winds and wind stress, the presence of 
ice cover, a small barotropic tide, and a weak cir- 
culation on basin scales, Levine et al. [1987]. Peaks 
at the M2 tidal frequency and its harmonics are not 
apparent. There was no discernible kinetic energy 
dissipation rate associated with the wavefield. Pad- 
man and Dillon [1987], consistent with low levels 
of energy and shear. 

The frequency spectra for both temperature and 
velocity were much whiter than typical Levine et 
al. [1987]. For frequencies much larger than /, the 
power law a~^ tends to r = 1.2, Levine [1990]. 

Velocity profiles were obtained with Expendable 
Current Profilers (XCPs). These data suggest a 
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Figure 30. Ocean Storms frequency spectra (blue and 
green lines) of horizontal kinetic energy and potential en- 
ergy from the main thermocline (200 m). Black curves 
represent fits of (21) with r — 2.20. The thick vertical 
line represents the buoyancy frequency cut-off. 
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vertical wavenumber bandwidth about ten times 
larger than in the GM model, D'Asaro and More- 
head [1991]. 

Several notes of caution are in order. First, the 
data were not taken as a random sample relative to 
the eddy field. The Canadian Basin is populated 
by relatively intense submesoscale vortices, and 
the AIWEX XCP survey was intentionally taken 
between such features. It is not clear whether 
this leaves the data set prone to straining effects 
as are evident in the MODE data set. Secondly, 
D'Asaro and Morehead [1991] argue for the pres- 
ence of quasi-permanent velocity finestructure at 
40 m vertical wavelength with relatively small lev- 
els of quasi-permanent density finestructure. Their 
diagnostics, though, depend upon being able to 
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Figure 31. Strex vertical wavenumber kinetic Ek- The 
overlay represents a fit of l/(mj + m^)^'^^ to the pre- 
frontal spectra, with m, = 0.0070 cpm (equivalent to 
j* = 3). Figure from D'Asaro [1984], their Fig. 4. Per- 
mission from J. Phys. Oceanogr., , copyright American 
Meteorological Society 1986. . 



account for all noise in the XCP measurements. 
Their noise budget does not include noise induced 
by fall rate variations of the XCP [Sanford et al., 
1993], which shows up as noise in the north-south 
velocity component at 15-30 m vertical wavelength. 
Finally, the spectral levels are so small that one 
might call into question whether nonlinearity is im- 
portant in shaping the spectrum. Boundary dissi- 
pation has been cited as a dominant decay process 
rather than interior wave breaking, Pinkel [2005]. 

4. SUMMARY 

Despite the fact that major deviations in the 
model parameters have been noted near boundaries 
[Wunsch and Webb, 1979] and near the equator 
[Eriksen, 1985a], it is astonishing that the char- 
acteristic shape, and after buoyancy scaling, level, 
are noted to be nearly universal in the literature. 
Variability in the spectral characteristics of the 
deep ocean internal wavefield is apparent. That 
variability, though, tends to be subtle. 

A number of patterns emerge from our analysis: 
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Figure 32. Alwex frequency spectra of isopycnal dis- 
placement gradients. The dashed line represents the GM 
frequency spectrum, E{a) ex a~^{a^ —f^)~^''^. The thick 
line represents E{a) ex a~'''^{a'^ — /^)~^" with with two 
values for /. (Figure extracted from Levine et al. [1987], 
their Fig. 2, permission from J. Geophys. Res., Copy- 
right 1987 American Geophysical Union). 
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• Geography: The variabihty of spectral shape 
is geographically related, rather than a function 
of spectral level. Measurements from Site-D ex- 
hibit characteristics similar to those from the east- 
ern Subtropical North Pacific and have power laws 
in good agreement with the GM76 model. Both 
are eddy-rich regions. Observations from the east- 
ern North Atlantic and the Arctic differ most sig- 
nificantly from the CM model. These regions 
have anomalously weak eddy fields and the inter- 
nal wave spectral levels are lower. Spectral shapes 
from the Sargasso Sea and the eastern boundary 
of the Subpolar North Pacific are intermediate in 
character. 

• Seasonality: A seasonal cycle, most promi- 
nent in high frequency energy, is apparent in the 
western Subtropical North Atlantic. A seasonal cy- 
cle is also apparent in the eastern Subpolar North 
Pacific. Their character, though, differs. The 
high frequency signal in the North Pacific lags 
the near-inertial cycle, while near-inertial and high 
frequency seasonal variability are in phase in the 
North Atlantic. The eddy fields in the two regions 
are distinctly different. 
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Figure 33. Alwex vertical wavenumber spectra of hor- 
izontal velocity. Thin lines represent fits of the GM75 
spectrum, Ek oc (f/m, + m)^'^, with different values of 
j, . This spectrum does not approach its asymptotic roll- 
off quickly. The thick overlay is Ek oc (l/{ml + m^))"-''^ 
with j, = 50 and total energy 0.95 x fO"'* m^ s"^. (Fig- 
ure extracted from D'Asaro and Morehead [f99f], their 
Fig. 7, permission from J. Geophys. Res., Copyright 
199f America Geophysical Union). 



• Amplitude: The GM spectral model has an 
energy of Eq = 30 cm^ s~^. With the stipulated 
frequency distribution (21), this tends to overes- 
timate observed frequency spectra from the deep 
ocean for a ^ f. This pattern is evident in tables 
presented by Wunsch [1976]; Wunsch and Webb 
[1979]; Fu [1981]; Nowlin et al. [1986]. Despite this 
overestimate, the total energy in the internal wave 
band agrees much better with Eq. The difference is 
additional energy in the tides and at near-inertial 
frequencies. 

• Separability: The observed deep ocean inter- 
nal wavefield is not, in general, separable. When 
available, vertical wavenumber spectra of poten- 
tial energy tend to have larger characteristic ver- 
tical scales than vertical wavenumber spectra of 
kinetic energy. Fits to the kinetic energy spectra 
produce roll-offs ranging from 3 < j* < 15 — 20, 
with lower amplitude frequency spectra being 
characterized by larger values of j*. Fits to the 
potential energy spectra produce j=k = 1 — 2. Non- 
separability is in the sense of a high-wavenumber 
near-inertial addition to a separable spectrum. In- 
dications of this are also apparent in the IWEX 
spectrum. The issue of separability appears to 
have been considered only in the continuum limit 
(/ < a < iV) by Garrett and Munk [1972]. 

• Covariability: Frequency spectra with power 
laws whiter than the GM model (o""^) tend to have 
vertical wavenumber power laws that are redder 
than the GM model {m~'^). 

• The mid- frequency dip. Many of the fre- 
quency spectra show a departure from the smooth 
spectral models at frequencies between 1.5 / and 
4-5 cph (e.g.. Fig. 14). Apart from tides and 
tidal/inertial peak harmonics, the observed spec- 
tra tend to be at somewhat lower levels than the 
model fits in this frequency band. Levine [2002] 
notes a similar pattern and suggests it is linked 
to internal wave dynamics at frequencies less than 
semi-diurnal. 

• Bandwidth: A recurring pattern is one of 
the near-inertial field having a larger vertical 
wavenumber bandwidth than the rest of the inter- 
nal wavefield. Data sets were selected on the ba- 
sis that they represent open ocean conditions and 
are not in the near-field of sites associated with 
topographic generation processes. We have con- 
sequently sought to relate the near-inertial band- 
width to interior processes. An unsettling possi- 
bility is that topographic infiuences may not be 
appreciated. 

The GM76 model is a reasonable description of 
the winter time internal wave spectrum at Site-D. 
The representativeness of the GM model for the 
background state here is not surprising as the GM 
frequency spectrum is no more and no less than a 
curve fit to data from Site-D. But what may not 
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be appreciated is that those original data were ob- 
tained during the dead of winter. At other times 
and in other places, differences between this sim- 
plistic model and observations are quite apparent. 
Departures from the GM model represent clues 
about the relative importance of the various possi- 
ble sources and how those source inputs are modi- 
fied by dissipation, wave/wave and wave-mean in- 
teractions. How those clues can be brought to- 
gether is discussed in the following section. 

5. SPECTRAL BALANCES AND RADIATION 
BALANCE EQUATIONS 

In the previous section we documented variabil- 
ity in the background internal wave spectrum. Our 
hypothesis is that such background states are in 
near equilibrium between forcing, dissipation and 
energy transfers. In this section we review what we 
perceive to be the key elements underlying such a 
hypothesis. We do not attempt an exhaustive re- 
view nor do we intend to elevate any result beyond 
the status of plausibility. As a prelude to that dis- 
cussion, we present a 35 year old paradigm for the 
spectral balances as a strawman. What we will find 
in the ensuing discussion is that key components 
of the strawman require reconsideration. 

Figure 34 from McComas [1977] presents one 
version of the internal wave energy balance in the 
spectral domain. A similar diagram can be in- 
ferred from Miiller and Olbers [1975]. Energy is 
gained in region I from surface waves at low verti- 
cal wavenumber and high frequency. In the over- 
lapping Region II (denoted S.D.), energy is input 
through interaction with the mesoscale eddy field 
and moved to region III. Nonlinear interactions 
transfer energy from Region III to Region IV via 
the Parametric Subharmonic Instability. Region V 
is in equilibrium with respect to nonlinear transfers 
associated with the Induced Diffusion mechanism 
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Figure 34. The McComas [1977] model for energy flow 
in the deep sea internal wavefield. 



and transports energy to Region VI, where dissi- 
pation occurs. 

We use this strawman to draw focus onto the 
following questions: 

• do the observed spectra represent approxi- 
mate steady states? 

• are the sources appropriately aligned along 
or within the spectral domain and of sufficient 
strength as to balance nonlinear transfers? 

• does variability of the spectral parameters re- 
late to variability in the forcing? 

Intrinsic difficulties with the McComas [1977] di- 
agram are that the source distributions are likely 
to be significantly different than stated and there 
are significant questions about the characterization 
of nonlinear transports associated with the weakly 
nonlinear theory that provides the named mecha- 
nisms above. 

We proceed by noting that a quantitative ex- 
planation of the observed spectral variability re- 
quires a corresponding analysis in the spectral do- 
main. In the case of media with weak inhomo- 
geneities the spectral evolution is described by a 
kinetic equation that has both transfer terms as- 
sociated with nonlinearity and a transport term 
that describes the motion of wave packets through 
the momentum-coordinate space. The transport 
term is a wave analog of Liouville's theorem or the 
continuity equations for a distribution function in 
statistical mechanics [Landau and Lifshitz, 1980; 
Pitaevskii and Lifshitz, 1981]. The kinetic equa- 
tion with the transport term is called a radiative 
balance equation, and is given by Lvov et al. [2003]: 



dn{p) 
dt 



+ Vr • [(u + Cg)n(p)] + Vp•[7^n(p)] 



Soip) - Si{p) +Tr{p). 



(23) 



The factor Cg = VpCj represents the group veloc- 
ity, u subinertial currents, 7^ = — Vr(i^ + p • u) 
refractive effects associated with spatially inhomo- 
geneous stratification and subinertial currents, T^ 
transfers of action, So interior sources and Si sinks 
(dissipation). Gradient operators in the spatial 
and spectral domains are Vr and Vp, respectively. 
Our underlying hypothesis is that the observed 
variability can be understood as spatially local sta- 
tionary states of the radiation balance equation 
(23). Thus we consider representations of nonlin- 
ear transfers Tr{p) in Section 5.1, a source func- 
tion So{p) characterization of wave-mean interac- 
tions {corresponding to terms Vr- [(u + Cg)n(p)] + 
Vp•[7^?^(p)]} in Section 5.2, external forcing So{p) 
in Section 5.3 and dissipation Si{p) in Section 
5.4. Possible issues related to propagation are ad- 
dressed in Section 5.5. See Fig. 35 for a schematic 
table of contents. 
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Wave Mean Interactions 

Section S.2 Eqn 70-73 Figs 44,45 



Surface-Internal Wave 
Coupling Section 5.3.2 



Tides 
Section 5.3.4 



Variable Wind Stress 
Section 5.3.1 



Coupling to Mixed Layer 
Section 53.3 



Wave Wave Interactions 

Section 5.1 Equation 37 and 59-«l Figures 37, 38, 40 

Approximations: eqns 41-42, 50-57 



Dissipation section 54 



Figure 35. A schematic Table of Contents for Section 5. The radiative balance equation describes the 
transfer energy in the spectral and spatial domain. Major processes are identified and classified with 
section, equations and figure numbers. 



5.1. Nonlinearity Trip) in (23) 

We began this investigation with the hypothesis 
that the observations, particularly the signature of 
covariable power laws, represent the end product 
of a nonlinear equilibration process. Nonlinear in- 
teractions between internal waves is a fascinating 
subject with a significant literature spanning 50 
years. An inclusive and insightful review is given 
in Miiller et al. [1986]. Our discussion below draws 
from this and includes more recent work. 

Our starting place is weak nonlinearity as de- 
scribed by a kinetic equation. A kinetic equation 
is a closed equation for the time evolution of the 
wave action spectrum in a system of weakly inter- 
acting waves. It is usually derived as a central 
result of wave turbulence theory. The concepts 
of wave turbulence theory provide a fairly general 
framework for studying the statistical steady states 
in a large class of weakly interacting and weakly 
nonlinear many-body or many-wave systems. In 
its essence, classical wave turbulence theory [Za- 
kharov et al., 1992] is a perturbation expansion in 
the amplitude of the nonlinearity, supposing linear 
plane wave solutions of ae*^'"'^"'^*^ at leading or- 
der and a slow amplitude modulation at the next 
order of the expansion, a(r)e*('"'P~'^*), by resonant 
interactions. This modulation leads to a redistri- 
bution of the spectral energy density among space- 
and time-scales. Below we sketch the derivation of 
such an equation following Zakharov et al. [1992]. 

Nonlinearities in the equations of motion are 
quadratic, so assuming the nonlinearity to be weak, 
first order expressions in those equations contain 
product terms of two waves having wavenumber- 
frequency (pi,(Ti) and (p2,o"2). For example, 



— Jl{pi+P2)-r-{ai+cr2)t] 



(24) 



If the combined wavenumber pi+P2 and frequency 
o"! + o"2 match those associated with a free wave, 
resonance occurs and the amplitude of the third 
wave grows linearly in time as long as its ampli- 
tude remains small. In general, the three wave 



resonance conditions can be written: 



p = Pi ± p2 and (T = 0"! ± (72 



(25) 



Solutions to (25) are referred to as the resonant 
manifold and the three waves form a resonant 
triad. The possibility of solutions to (25) depends 
simply upon the geometry of the dispersion sur- 
faces [Phillips, I960]. Such three wave solutions 
are possible for internal waves but not, for exam- 
ple, surface gravity waves in deep water. 

The evolution of wave amplitude ap follows most 
directly from Hamilton's equation: 



.d 

dt ^ 



6n 



(26) 



with Hamiltonian Ti that is nominally the sum of 
kinetic and potential energies, or Lagrange's equa- 
tion: 



d dC dC 



dt dcbr 



duv 







(27) 



with Lagrangian C that is the difference of kinetic 
and potential energies. 

Usually, but not always, one needs to adopt a 
linearization to obtain the Hamiltonian Ti in (26) 
or the Lagrangian C in (27). The difficulty is that, 
in order to utilize (26) or (27), the Hamiltonian or 
Lagrangian must first be constructed as a func- 
tion of the generalized coordinates and momenta 
(op). It is not always possible to to do so explic- 
itly, in which case one must set up a perturbation 
expansion which imposes a small amplitude limita- 
tion upon the result. See Section 5.1.3.1 for further 
detail. 

The Hamiltonian of a system with quadratic 
nonlinearity can be expressed as [Zakharov et al., 
1992]: 



n 



dp a 



Pl^pl 
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+ '^-p + Pl+P2(^p^l,p2«p«Pl«P2 +C-C.) 



(28) 



after a Fourier decomposition. The interaction co- 
efficients U and V are extended algebraic expres- 
sions involving wave amplitude Op, frequency dp 
and momentum p of the three waves. 

Having obtained the Hamiltonian (28) one intro- 
duces wave action as 



rir 



{a a 



p^p/) 



(29) 



where (...) means the averaging over statistical en- 
semble of many realizations of the internal waves. 
To derive the time evolution of n^ the amplitude 
equation (26) with Hamiltonian (28) is multiplied 
by a*, the amplitude evolution equation for a* is 
multiplied by a, the two equations are differenced 
and the result averaged (. . .) to obtain: 



9np 

dt 



^P^1P2'^^1P2'^(P-P1-P2) 



-^pPpVpp/(P2-P-Pl) 
-^P%^%'^(P1-P2-P) 



(ipi(ip2, (30) 



where we introduced a triple correlation function 



•^PiP2'^(Pi - P - P2) = (apapiap2). 



(31) 



If we were to have non-interacting fields, i.e., fields 
with V^i-p2 being zero, this triple correlation func- 
tion would be zero. We therefore invoke a per- 
turbation expansion in smallness of interactions to 
calculate the triple correlation at first order. The 
first order expression for dn-pjdt thus requires com- 
puting dJpip2/dt to first order. To do so we take 
definition (31) and use (26) with Hamiltonian (28) 
and apply (...) averaging. We get 
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P3V 



/PI 
P2P3 



= / [-^(Vp%s)*«l'^(Pl - P4 - P5) 

+(^p^.P5r«:^(p4-P2-P5) 

+^p''3P5'^P2P4''^(P4 - P3 - P5)] dpidli^2) 

Here we introduced the quadruple correlation func- 
tion 

■^Pspf '^(Pi + P2 - P3 - P4) = {aiyp^ap-^ap^) .(33) 

The next step is to assume Gaussian statistics, and 
to express Jp3P4 as a product of two two-point cor- 
relators as 



Then 
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P2 
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P3^ 



TPl 
P2P3 



P2P3^ 



(nins + nin2 - n2n3) 



(34) 



Time integration of the equation for Jp2P3 will 
contain fast oscillations due to the initial value 
of Jp2p3 and slow evolution due to the nonlinear 
wave interactions. Contributions from the first 
term will rapidly decrease with time, so neglect- 
ing these terms we get 
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+ ir. 



(35) 



'PI "^P2 "P3 ^ '-^ P1P2P3 

Here we introduced the nonlinear damping of the 
waves rp^p2p3 that relates to the breadth of the 
resonant manifold. We will elaborate on rpjp2P3 
in Section 5.1.3. We now substitute (35) into (30), 
assume for now that the damping of the wave is 
small, and use 



We 
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then obtain 
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-7r(5(A). 



(36) 



.A + ir_ 

the three-wave kinetic equa- 
tion [Zakharov et al., 1992; Lvov and Nazarenko, 
2004; Lvov et ai, 1997]: 
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Pl"P2 



'^(Pl-P3)'5(P2-P4)+<J(Pl-P4)(5(P2 



7pl2'5p-pi-p2'5(o-p - CTpi - 0-p2)dpi2 
^P2'pl^/l2p'^Pi-P2-p'^(c^Pi - ^P2 - 0-p)(ipi2 

-4vr/|FP^jV2pi Sp.,-p-p-,5{ap2 - fJp - ap^)dpi2 , 

with /pi2 = np^np2 - np(np^ + "'P2) • (37) 

The 6 functions require the wavenumbers and fre- 
quencies sum to zero and ensure that pseudomo- 
mentum [pn(p)] and energy [(Tpn(p)] spectral den- 
sities are conserved. The Resonant Interaction Ap- 
proximation (RIA) is represented by the reduction 
of (35) to the resonant manifold (25) and S function 
representation in (36). 

The typical assumptions needed for the deriva- 
tion of kinetic equations are: 

• Weak nonlinearity, 

• Gaussian statistics of the interacting wave 
field in wavenumber space and 

• Resonant wave-wave interactions 

In a more systematic derivation, such as Lvov and 
Nazarenko [2004], the kinetic equation can be ob- 
tained using only an assumption of weak nonlin- 
earity. A theory of weak interaction has its limita- 
tions and this motivates discussions of renormaliza- 
tion attempts to account for increasing nonlinear- 
ity (-,5.1.3.2) and ray-tracing descriptions of wave- 
~\fwe -interactions (Section 5.1.3.3). 
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5.1.1. Resonant Interactions: Scale Invari- 
ant Solutions 

Following Kolmogorov's viewpoint on energy 
cascades in isotropic Navier-Stokes turbulence, one 
may look for statistically stationary states us- 
ing scale-invariant solutions^ to the kinetic equa- 
tion (37). The solution may occur in an iner- 
tial subrange of wavenumbers and frequencies that 
are far from those where forcing and dissipation 
act, and also far from characteristic scales of the 
system, including the Coriolis frequency resulting 
from the rotation of the Earth, the buoyancy fre- 
quency due to stratification and the ocean depth. 
Under these assumptions, the dispersion relation 
and the interaction matrix elements are locally 
scale-invariant and solutions of (37) take the form 



n(k, m) 



\m\ 



(38) 



Values of a and b such that the right-hand side of 
(37) vanishes identically correspond to steady solu- 
tions of the kinetic equation, and hopefully also to 
statistically steady states of the ocean's wave field. 
To connect with the observations, the corre- 
sponding vertical wavenumber - frequency spec- 
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Figure 36. The observational points. The filled cir- 
cles represent the Pelinovsky-Raevsky (PR) spectrum 
[{a,b) = (3.5,0.5)], the convergent numerical solution 
[{a,b) = (3.7,0.0)] determined in Lvov et al. [2010] and 
the CM spectrum [{a,b) = (4.0,0.0)]. Circles with stars 
represent estimates based upon one-dimensional spec- 
tra from the western North Atlantic south of the Gulf 
Stream (IWEX, FASINEX and SFTRE/PMIII), the east- 
ern North Pacific (STREX/OS and PATCHEX^), the 
western North Atlantic north of the Gulf Stream (Site- 
D), the Arctic (AIWEX) and the eastern North Atlantic 
(NATRE^ and NATRE^). There are two estimates ob- 
tained from two-dimensional data sets from the eastern 
North Pacific (SWAPP and PATCHEX^) represented as 
circles with cross hairs. NATRE^ and NATRE^ represent 
fits to the observed spectra and observed minus vorti- 
cal mode spectra, respectively. Lines of constant vertical 
wavenumber power laws for the observed energy spec- 
trum run diagonally downward. Lines of constant fre- 
quency power laws are vertical. 



trum of energy is obtained by transforming nk,m 
from the three dimensional wavenumber space 
(k, m) to the two dimensional vertical wave- 
number-frequency space ((T, m) and multiplying by 
frequency. In the high frequency-high wavenumber 
limit, 



E{7n, a) 



oc m 



-b^2- 



The frequency-wavenumber power laws inferred 
from the observations (m^'^o"^'') are thus related 
as: 

r = a — 2, q = a + b — 2. 

Parametric fits to the observations (Fig. 36) pro- 
duce a clustering about two regions of parameter 
space with several notable outliers. Observations 
from Site-D and from the subtropical North Pa- 
cific tend to cluster about the GM76 power laws 
(m~^o"~^). These data were obtained in regions 
poleward of the critical latitude (28.9) at which the 
M2 semidiurnal tide can decay into near-inertial 
frequencies through a parametric subharmonic in- 
stability. The IVIATE data from the subpolar North 
Pacific and observations from south of the Gulf 
Stream are grouped together with power laws de- 
scribing spectra distinctly redder (whiter) in the 
vertical wavenumber (frequency) domain. This 
grouping includes data both poleward and equa- 
torward of the semidiurnal critical latitude. Obser- 
vations from the eastern North Atlantic (NATRE) , 
the Arctic (AIWEX) and, tentatively, the subpolar 
North Pacific (Ocean Storms) appear as outliers. 
5.1.1.1. Nonrotating (/ = O) limit 

Unlike Kolmogorov turbulence, the exponents 
which give steady solutions to (37) can not be de- 
termined by dimensional analysis alone. This is the 
case owing to multiple characteristic length scales 
in anisotropic systems, e.g., [Polzin, 2004a]. The 
anisotropic analysis is much more complicated and 
as a prelude to that discussion we enquire whether 
the improper integrals in (37) converge. This is re- 
lated to the question of locality of the interactions: 
a convergent integral characterizes the physical sce- 
nario where interactions of neighboring wavenum- 
bers dominate the evolution of the wave spectrum, 
while a divergent one implies that distant, nonlocal 
interactions in the wavenumber space dominate. 

It turns out that the internal-wave collision in- 
tegral diverges^ for almost all values of a and b. 
In particular, the collision integral has an infra-red 
(IR) divergence at zero, i.e., |ki| or |k2| — )• and 
an ultra-violet (UV) divergence at infinity, i.e., |ki| 
and |k2| — )• 00. A detailed analyses performed in 
[Lvov et al., 2010] returns the convergence condi- 
tions in the IR limit, ki or k2 — )• of: 



a -h 6/2 - 7/2 < and - 3 < 6 < 3, 
a - 4 < and 6 = 0, 
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or 



a - 7/2 < and b= 1, 
a + 6 — 5 < and 6 > 3, 

a — 5 < and b < —3. 



(39) 



Similarly, convergence in the UV limit 
and I k2 |— )■ oo implies that 

a + 6/2 - 4 > and - 2 < 6 < 2, 

a - 7/2 > and b = 0, 

a - 3 > and 6 > 2, 



or 



a + 6 - 3 > and b < -2. 



(40) 



The domains of divergence and convergence are 
shown in Fig. 37. 

There is only one exception where the integral 
converges in both limits: the line segment along 
6 = with 7/2 < a < 4. Moreover, there is one 
special value, (a, 6) = (3.7,0), in which IR and 
UV contributions are both convergent and cancel. 
That is, (a, 6) = (3.7,0) represents a stationary 
state state of the kinetic equation in the absence 
of rotation. Note that the line 6 = is 'special' in 
that it corresponds to wave action independent of 
vertical wavenumbers, dn{p)/dm = 0. This con- 
vergent solution lies in the general proximity of 
GM76 (a, 6) = (4, 0) and observations from Site- 
D, SWAPP and PATCHEX. 

The possibility also exists that one can construct 
a divergent solution in the following sense. If 
the UV and IR nonintegrable singularities have 
opposite signs for some (a, 6), it may be possi- 
ble to devise a conformal mapping that matches 
the oppositely signed singularities. In this case 
one would have, at least on the level of a techni- 
cal mathematical exercise, an approximate station- 
ary state. Pelinovsky and Raevsky [1977]; Zeitlin 
[1992]; Lvov and Tabak [2001] find this to be pos- 
sible for (a, 6) = (3.5,0.5). Lvov et al. [2010] de- 
termine that such solutions are more general, and 
possible approximate stationary states occupy the 
dark grey shaded regions of Figure 37. The black 
regions are not permitted due to UV and IR diver- 
gences having similar signs. 

The parametric spectral estimates tend to have 
larger a values than either the convergent or the 
Pelinovsky and Raevsky [1977] solution and also 
tend to lie on the border of UV convergence. One 
instance (NATRE) is seen to occupy the black re- 
gion, but this point migrates to a state of UV con- 
vergence if an attempt is made to subtract non- 
wave finestructure (See Section 3.3.4 for details). 

Although it represents a possible stationary 
state, a balance between divergent and oppositely 
signed integrals is not a satisfactory physical sit- 



uation. What it represents is a statement that 
integration endpoints other than ±00 need to be 
inserted in (37). Obvious candidates are / in the 
IR limit and A^ in the UV limit. 

5.1.1.2. / ^ 0: Extreme Scale Separated In- 
teractions 

When the kinetic equation is evaluated on the 
resonance surface, extreme scale separated interac- 
tions are believed to dominate the transfers [Lvov 
et al., 2010; Miiller et al., 1986]. Thus approaching 
the issue of scale invariant solutions to the kinetic 
equation is not as simple as a spectrally local en- 
ergy cascade. A detailed analysis of the extreme 
scale separated interactions is required. Three sim- 
ple interaction mechanisms were identified by Mc- 
Comas and Bretherton [1977] in the limit of an 
extreme scale separation. The limiting cases are: 

• the vertical backscattering of a high-frequency 
wave by a low frequency wave of twice the vertical 
wavenumber into a second high-frequency wave of 
oppositely signed vertical wavenumber. This type 
of interaction is called elastic scattering (ES). 

• The scattering of a high-frequency wave by 
a low-frequency, small-wavenumber wave into a 
second, nearly identical, high-frequency large- 
wavenumber wave. This type of interaction is 
called induced diffusion (ID). 

• The decay of a low wavenumber wave into 
two high vertical wavenumber waves of approxi- 
mately one-half the frequency. This is referred to 
as the parametric subharmonic instability (PSI). 
This mechanism was identified by McComas and 
Bretherton [1977] as causing the transfer of energy 
from frequencies of 2/ and low vertical wavenum- 
ber to high vertical wavenumber near inertial os- 
cillations in the GM spectra. 

Below we review how these extreme scale sepa- 
rated interactions can address the observed covari- 
ability in power law fits. 

5.1.1.3. Induced Diffusion and covariable 
po^ver laws 

One can further reduce (37) in the ID limit to a 
Fokker-Plank equation [Gardiner, 2004; McComas 
and Bretherton, 1977] 



dn{p) 
dt 



dpi dpj 



(41) 



such that 



n 



V 



d^p'd^p" 



n(p")p'/p-|^p^,p" 



'6{p-p'-p")5{a-a'-a'm 



in which p and p' are the wavevectors of high fre- 
quency waves, p" is the wavevector of the low fre- 
quency member of the triad and \V^, //[^ = ^a"k'^. 
To obtain the Fokker-Plank equation one: (i) ne- 
glects terms quadratic in high frequency amplitude 
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relative to linear terms, (ii) expands [n(p — p") — 
n(p)] in a Taylor series, (iii) sums the formulae for 
the two cases | p' |>| p | and | p' |<| p |, and, not- 
ing the subtractive cancelation between the two, 
performs a second Taylor series expansion. McCo- 
mas and Bretherton [1977] present scale estimates 
of (41) and argue that the vertical coordinate dom- 
inates the transfers, thereby reducing the integral 
expression to being proportional to the vertical 
gradient variance (i.e., shear variance) along the 
resonance curve. We reconsider these arguments 
in Section 5.1. I.e. 

Substituting the scale invariant solution (38) 
into (41), one obtains, for the GM class of spec- 
tra in which near-inertial frequencies dominate the 
shear variance, stationary states if either 



9 - 2a - 36 = or 6 = 0, 



(43) 



Lvov et al. [2010] obtain this result via rigorous 
asymptotic arguments, independent of the reduc- 
tion of the kinetic equation (37) to the Fokker- 
Plank equation (41). McComas and Miiller [1981b] 
interpret 6 = as a no action flux in vertical 
wavenumber domain, while 9 — 2a — 36 = is 
a constant action flux solution. Power laws 6^0 
imply a down-gradient action flux to higher/lower 
vertical wavenumber along lines of constant hor- 
izontal wavenumber. Associated with this is the 
transport of action to lower/higher frequency and 
the gain/loss of energy by the near-inertial field. 
A positive residual to (37) implies an increase in 
high vertical wavenumber, high frequency spectral 
density and is consonant with a balance between 
resonant transfers and a high wavenumber sink. 

The situation changes in the absence of rotation 
(/ = 0). In this case, similar manipulations sug- 
gest that stationary states can be obtained if either 



6-1 = or 6 = 0. 



(44) 



in which 6 = 1 is the constant action flux solution 
of the Fokker-Plank equation. 

The majority of the observational data points 
cluster in two groups that are not far from the ID 
stationary states for scale invariant solutions (Fig. 
37). The observations exhibit a tendency to lie in 
regions having either small transport (6 = 0) or 
transport to higher vertical wavenumber and lower 
frequency (6 > 0). The observations also tend to 
lie in regions having a positive residual, for which 
^ > 0, and is consistent with a balance between 
the ID flux convergence and an unspecified sink. 
The most extreme departure from the ID station- 
ary states is NATRE, which is south of the critical 
latitude for the semi-diurnal lunar tide and has an 
anomalously large vertical wavenumber bandwidth 

ij* = 18). 

However appealing the apparent pattern match 
between the observations and stationary states un- 
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horizontal power law of 3-d action spectrjm (a) 



Figure 37. The observational points and the theories 
[Lvov et aL, 2010]. The filled circles represent scale- 
invariant stationary states identified by Pelinovsky and 
Raevsky [1977], a convergent numerical solution deter- 
mined in Lvov et al. [2010], the GM spectrum and two 
possible dynamic balances identified in McComas and 
Miiller [1981b]. Circles with stars represent power law 
estimates based upon one-dimensional spectra. Circles 
with cross hairs represent estimates based upon two- 
dimensional data sets. See Fig. 36 for the identifica- 
tion of the field programs. Light grey shading represents 
regions of the power-law domain for which the collision 
integral converges in either the infrared or ultraviolet 
limit. The dark grey shading represents the region of the 
power law domain for which neither the infrared nor 
ultraviolet limits converge. The region of black shading 
represents the subdomain for which the infrared and ul- 
traviolet divergences have the same sign. Additionally, 
a line segment between 3 < a < 7/2 along b = 1 is 
IR convergent and a line segment between 7/2 < o < 4 
along 6 = is IR and UV convergent. These line seg- 
ments arise as higher order cancellations in the integrand. 
Overlain as solid white lines are the induced diffusion sta- 
tionary states. Thin white contours are proportional to 
the factor 6(9 — 2a — 3b) in (41) with contour intervals of 
±(0.25,0.50,0.75,1.0,2.0,3.0,...). The +- symbols de- 
note the sign of the action tendency (37) in the a — b 
domain. 



der the ID mechanism may be, we believe there is 
more to the story. That story includes numerical 
evaluations of the kinetic equation which contain 
significant residuals inconsistent with G1VI76 being 
a stationary state. These residuals lead us to con- 
sider how increasing levels of nonlinearity at high 
wavenumber could change our interpretations in 
Section 5.1.3. Before doing so, though, we fill in 
a bit of the background by discussing the dynamic 
balance of IvIcComas and IVIiiller [1981b]. 
5.1.1.4. The dynamic balance 

A quantitative formulation of the internal wave 
energy budget is given structure by the dynamic 
balance put forward in IVIcComas and IVIiiller 
[1981b]. This key study ascribes an inertial sub- 
range character to the G1VI76 spectrum after inte- 
grating over the frequency domain. The dynamic 
balance is constructed by dividing the frequency 
domain into distinct regions and characterizing 
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each region with transfers associated with the dom- 
inant triad class. Thus frequencies of / < o" < 2/ 
and vertical wavenumber m gain energy via PSI 
transfers from frequencies of 2/ < o" < 4/ and 
lower wavenumber. Frequencies of 2 < a < 4/ and 
wavenumber m loose energy via PSI transfers to 
frequencies of / < o" < 2/ and higher wavenumber 
and also gain energy via ID energy spectral fluxes 
across a = 4/. Frequencies o" > 4/ and wavenum- 
ber m gain or loose energy as action is transfered 
to higher or lower vertical wavenumber in associ- 
ation with ID coupling to near-inertial waves of 
smaller vertical wavenumber. Finally, frequencies 
a = f and wavenumber m gain or loose energy 
in association with ID coupling to high frequency 
waves a > Af . These domains are apparent in the 
schematic from McComas [1977] (Figure 34) and 
in numerical evaluations of the kinetic equation, 
e.g.. Figure 40. Matching expressions for the en- 
ergy spectrum and its time rate of change across 
the boundary o" = 4/ provides the high frequency 
solution: 



Boltzman rate: 



Eid{(t,P) = E 



-4[1 + |^0 (45) 
a^ QAx 4/ 



/32 



which is nicely fit by a power law n(p) oc 
. Matching along /3 = /?* provides 



k 



-3.875^-0.125 



'm 



Eid{(t,P) 






64x'''/3j 



(46) 



-4.0^-0.20 



which is fit by a power law n(p) oc k 
Both solutions are plotted in Fig. 37. 

One of the attractive properties of the first solu- 
tion [(a, 6) = (3.875,0.125)] is that the expressions 
for energy transport to small vertical scale lead 
to estimates that are difficult to distinguish from 
turbulent dissipation estimates obtained as part of 
the PATCHEX field program [Gregg and Sanford, 
1988]. These estimates, along with those obtained 
in NATRE [Polzin et al., 1995], pin down the GM76 
dissipation rate at approximately 8 x 10"^" W/kg. 
See Polzin [2004a] for discussion. 

There are, however, several disquieting issues. 
The first is that the downscale energy transport 
occurs only in vertical wavenumber and thus re- 
quires an energy source at high frequency. The 
second is the requirement of integrating over the 
frequency domain and invocation of dissipation to 
obtain a stationary state. The RIA dissipation 
mechanism (Section 5.4.2) acts primarily on near- 
inertial shear whereas numerical evaluations of the 
kinetic equation (Section 5.1.1.e) indicate signifi- 
cant non-stationary tendencies at high frequencies. 
5.1.1.5. Numerical Evaluations of (37) 

The analysis presented above was based upon 
analytic work. Here we report numerical evalu- 
ations of the kinetic equation expressed via the 



l^B 



HP) 
2n(p) 



(47) 



The Boltzman rate represents the net rate of trans- 
fer for wavenumber p and is a low order measure 
of nonlinearity for smooth, isotropic and homoge- 
neous spectra. The individual rates of transfer into 
and out of p maybe significantly larger for spec- 
tral spikes [Pomphrey et al., 1980; Miiller et al., 
1986] and potentially for smooth, homogeneous but 
anisotropic spectra. 

A characterization of the level of nonlinearity 
is made by normalizing the Boltzman rate by the 
characteristic linear time scale Tp = 27r/o"p: 



2Trhr 



(48) 



A stationary state is defined by ep = 0. The nor- 
malized Boltzman rate (48) serves as a low order 
consistency check for the various kinetic equation 
derivations. An 0(1) value of ep implies that the 
derivation of the kinetic equation is internally in- 
consistent. 

We present evaluations of the kinetic equation 
as discussed in Lvov et al. [2011]: we use a numer- 
ical scheme designed for off resonant calculations 
using the Lvov and Tabak [2004] Hamiltonian and 
then we take the resonant limit. See also Section 
5.1.3.1. We use two spectra that mimic those at 
Site-D [GM76 with (a, b) = (4.0, 0.0)] and the Sar- 
gasso Sea [(a, 6) = (3.75,0.50), Appendix D]. 

The characteristic pattern for GM76 is two 
positive and one negative lobes in the vertical 
wavenumber - frequency domain with boundaries 
at approximately 2/ and 5/ (Fig. 40). Action loss 
at frequencies 2/ < o" < 5/ is qualitatively con- 
sistent with the PSI mechanism transferring en- 
ergy into near-inertial frequencies, [McComas and 
Miiller, 1981a]. The residual at o" = 3/ is, how- 
ever, nearly an order of magnitude larger than 
what would be predicted from the PSI time scale 
[McComas and Miiller, 1981a], which leads to: 



PSI 



8-ir^ 27 f m 
lOx 32 a rrir 



(49) 



with X = vlO and rric a high wavenumber cut- 
off defined in (17). Comparative figures for GM76 
can be found in McComas and Bretherton [1977] 
(Fig. 10), McComas [1977] (Fig. 3), McComas 
and Miiller [1981b] (Figure 11) and Pomphrey et 
al. [1980] (Fig. 6). 

Normalized Boltzman rates within the ID regime 
(T > 5/ are both 0(1) and positive. The fact that 
they are 0(1) is indicative that GM76 is not a 
stationary state, and the fact that they are posi- 
tive is troubling in that the nearby state of (45) 
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was specifically selected by McConias and Miiller 
[1981b] as having ID transfers with n(p) < to 
match a similar signed tendency along the ID/PSI 
boundary a = 4/. To resolve this we interpret the 
results in terms of the convergence conditions of 
the collision integral. 

The rigorous asymptotic analysis presented in 
Lvov et al. [2010] finds that both ultraviolet and 
infrared limits to the collision integral converge 
along 3.5 < a < 4 and 6 = 0. Thus interactions 
along this line segment are local in the nonrotating 
limit. With the addition of rotation the mathe- 
matical properties of the collision integral could be 
altered, but note that the line 6 = corresponds to 
an action spectrum with vanishingly small ampli- 
tude gradients in vertical wavenumber, with corre- 
spondingly small transfers that depend upon those 
gradients, such as the ID mechanism. We therefore 
infer that the positive contours tending vertical at 
the highest frequencies are unrelated to extreme 
scale separated ID transfers. 

We double check this inference by noting that lo- 
cality implies frequencies of / < o" < 2/ make little 
contribution to the collision integral. If such low 
frequencies are excluded from the numerical evalu- 
ations, we find n(p) to be virtually unchanged for 
GM76 at frequencies a > 4/ (Fig. 38). Similarly, 
we find little change if the wavenumber domain is 
limited to m < 2 rric and k < 2 rric- The infiuence 
of extreme scale separated ID transfers in Fig. 38 
would be apparent as contours of Cp tending to 
be parallel to a/m = constant in the vicinity of 
m = m^a/f, i.e., at small vertical wavenumbers 
prior to the asymptotic roll-off. Such an infiuence 
is not apparent in any of these GM76 evaluations. 

Finally, we note that, in the asymptotic limit of 
high vertical wavenumber, the gradients of the 3-d 
action spectrum in vertical wavenumber are van- 
ishingly small for 6 = 0. We therefore hypothesize 
that local interactions giving rise to these residuals 
are related to the gradients of the 3-d action spec- 
trum in horizontal wavenumber. A consequence of 
transfers in horizontal wavenumber is that it poten- 
tially changes the energy source at high frequency 
required by the vertical wavenumber Fokker-Plank 
equation (41)! 

The Sargasso Sea power law combination repre- 
sents a constant action flux solution 9 — 2a — 36 = 
to the Fokker-Plank equation (41). Comparative 
figures for the nearby state of GM75 [(a, 6) = 
(4.0,0.5)] can be found in Olbers [1976] (Fig. 4), 
McComas and Bretherton [1977] (Fig. 12) and 
Pomphrey et al. [1980] (Fig. 6). A bandwidth 
of jjf = 10 has been chosen as being representative 
of the region and the amplitude has been speci- 
fied to return a gradient variance at A^ = 10 m 
similar to the GM76 spectrum. The pattern for 
the Sargasso Sea is, excepting a = N, a, consis- 



tent negative tendency for a > 2f and for con- 
tours to be parallel to the ID resonance condition 
of m* = fm/a. Tendencies lead to ep = 0(1) at 
m = rric- Such 0(1) tendencies are inconsistent 
with the power law combination being a station- 
ary state, but we note that this is likely an artifact 
of the high wavenumber, high frequency domain 
{m, a) = (mc, N) not being within the asymptotic 
regime of a/m » f /m^,. 

Truncating the spectral domain to exclude fre- 
quencies o" < 2/ changes the tendencies consider- 
ably, so that the sign and magnitude are similar 
to those of GM76. Truncating the spectral domain 
to exclude vertical wavenumbers ra > 2mc changes 
only the results at highest vertical wavenumber. 
These diagnostics suggest the Sargasso Sea re- 
sults fit much better into the extreme scale sep- 
arated scenario, but local interactions cannot be 
neglected. 

DELETE THE PARAGRAPH BELOW: AIWEX 
Note the proximity to the line segment 3 < a < 3.5 
along 6 = 1. This line segment is an IR convergent 
state in the nonrotating analysis. We have not in- 
vestigated interactions in detail for this spectrum. 
5.1.2. Resonant Interactions: Departures 
from Scale Invariance 
5.1.2.1. NATRE: PSI decay of a tidal peak 

The parametric subharmonic instability de- 
scribed in Section 5.1.1.b has also been cited as 
providing a pathway for the transfer of energy from 
low-mode internal tides [Hibiya et al., 2002; MacK- 
innon and Winters, 2005]. The question of whether 
such PSI transfers can explain the unique character 
of the NATRE spectrum is examined in Polzin et 
al. [2008]. On the assumptions that m! <^ m,m" 
and cj' ~ 2o" ~ 2a" ^ 2/, the transfer integral 
reduces to 

Trip) = 2/dp'|I^p^,(p,_p)|'5[a -a' + a(p' - p)] 

{n(pO[n(p) + n(p' - p)] - n(p)n(p' - p)} (50) 

for the gain at high wavenumber. The transfer 
function | V^, , , n | ^ in this limit is given by [Mc- 
Comas and Miiller, 1981a]: 
,/3 



p',(p'-p)l 8 a{a'-a) 



a 



/2 



f 



a 



a 



2 a 



/2 



7V2 



-m ' 



'2 



(51) 



McComas and Miiller [1981a] note that with 
n(p) ~ ^^(p") ^ ^(pO the transfer integral fur- 
ther reduces: 
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Figure 38. Normalized Boltzman rates (48) calculated via (59) using the Lvov and Tabak [2004] ma- 
trix elements with S = 10~^. Left panels use the GM76 (18) spectrum, right panels use Sargasso Sea 
power laws (Appendix D). Upper panels utilize a spectrum that extends to m = 8mc, k = SrUc, with 
nic defined in (17). Middle panels: as above, but with E{a < 2f,m) = 0.0. Lower panels: as in the 
upper panels, but with restricted vertical m < 2mc and horizontal k < 2mc wavenumbers. Results are 
only plotted for m < rric. A value of 5 = 10~^ represents the resonant interaction limit. The black 
diagonal lines represent the ID resonance curve (//m* — a/m) extending from m,. Contour values are 
(0, ±0.01, ±0.02, ±0.05, ±0.1, ±0.2, ±0.5, ±1.0). 



with n{a, m) a 2-D action spectrum and semidi- Comas and Miiller [1981a] assume n(p') to be a 

continuous spectrum and here n(p') represents an 
internal tide, bandwidth hmited to a narrow range 

McComas and Miiller [1981a] and (52) is that Mc- of frequencies. Assuming further that the internal 



urnal tidal frequency gt- The difference between 
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tide is concentrated in niode-l (mi), the character- 
istic growth rate at high wavenumbers is exphcitly 
given by 



Tr{p)/n{p) 
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with E{a) the 1-d energy spectrum. 

The total rate of energy transfer from the semi- 
diurnal tide to the near-inertial field is 



aTr-d^k = / d^kE{p)T 



[P) 



Em/T, (54) 



with near-inertial energy Eni and time scale r from 
(38). This assumes that the near-inertial energy 
Eni is resonant with the tide. One of the restric- 
tions of the kinetic equation is apparent in (38). 
The tide can not be regarded as a single plane 
wave. Doing so implies £'((Tt) = Eo6{a — ax), so 
that the transfer rate is infinite. Numerical simu- 
lations of PSI transfers [MacKinnon and Winters, 
2005] indicate large and perhaps unrealistic trans- 
fer rates if forced with a single plane wave. 

When used as a diagnostic, Polzin et al. [2008] 
report that (54) tends to overestimate observed 
dissipation rates. An additional concern is that 
the Natre spectrum is unique in both its vertical 
wavenumber domain power law and extreme band- 
width. Implied is a scale selection process that is 
not represented in (54), which states that produc- 
tion is proportional to the shear in the tide times 
the near-inertial energy. 

This can be rectified by noting that the limit 
n{p),n{p") <^ n{p') does not need to be taken. 
Since n{p") = n(— p — p'), if one simply assumes 
that n(p") = n(— p), then 
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in which 7 is the minimum ratio of m' /m permitted 
by resonance constraints. The factor 7 is approx- 
imately 1.5-2.0 for Natre. The term quadratic in 
near-inertial amplitude represents a scale selection 
principle. 

A diagnostic for application to observed fields 
is obtained by summing the wave-action tendency 
equations for n(p) and n(p"), converting from a 
wavevector representation of the 3-D wave-action 
spectrum to a 2-D energy spectrum in vertical 
wavenumber-frequency coordinates and remember- 



ing that the observations, in distinction to the the- 
ory, reference a one-sided, rather than two-sided, 
spectrum. The conversion from 3-D to 2-D as- 
sumes that the wavefield is horizontally isotropic, 
result is that: 
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The total transfer from the tide to the near- 
inertial field is [Polzin et al., 2008]: 
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(58) 

Polzin et al. [2008] report a close correspondence 
between production (58) and the observed dissipa- 
tion rate, and the roll-off of the vertical wavenum- 
ber spectrum as a balance between 



m'^Eiar) 



and 



m 



E(a, m) 



5.1.2.2. The Dip: ID decay of a tidal peak 

The same argument used to provide a scale se- 
lection principle in PSI can be invoked for ID. 
The term quadratic in high frequency amplitude 
is dropped in the Induced Diffusion approximation 
(41). This term represents a consistent energy loss 
from the high frequency field to the low frequency 
tide. In this instance the quadratic term is largest 
when the scale separation is smallest and the argu- 
ments used to produce the ID approximation break 
down. Furuichi et al. [2005] report an analysis of 
energy transfers between a tide and a background 
sea of internal waves initialized as the GM spec- 
trum. They find an energy loss from the tide to the 
high frequency field along the ID resonance curve 
at relatively high frequencies and an energy gain by 
the tide from the high frequency field along the ID 
resonance curve at relatively low frequencies (small 
scale separations) (Fig. 39). 

Such a pattern of gain and loss is consistent with 
the mid-frequency dip being a product of nonlin- 
ear interactions with the tide, as conjectured by 
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Figure 39. Nonlinear energy transfer rates from the lowest vertical wavenumber M2 internal tide, with 
blue shades representing energy transfer from the tide to the background wavefield [Furuichi et al., 2005] 
. The red circles depict the PSI (upper left) and ID (lower right) resonant curves. Diagonal lines are 
curves of constant wave period. 



Levine [2002]. The essential question is whether 
this energy loss along the ID resonance curve is 
sufficiently large in comparison to other effects. 
5.1.2.3. The Dip: Departure from scale in- 
variance 

If the pattern of loss and gain in Fig. 38 is time 
stepped in the radiation balance equation, it will, 
at least initially, tend to result in a pattern resem- 
bling the mid-frequency dip, e.g.. Fig. 14. Much 
depends upon the magnitude and direction of non- 
linear transports and distribution of sources within 
the spectral domain / along the spectral bound- 
aries. It would help to understand how the mid- 
frequency dip is geographically distributed and re- 
lated to power law combinations in order to under- 
stand whether the mid-frequency dip results from 
nonlinearity. 

5.1.3. Resonance Broadening and Coordi- 
nate Representations 

Numerical evaluations of (37), Fig. 38, suggest 
that neither GM76 nor the Sargasso Sea spectrum 
is a stationary state, and this poses real problems 
for our hypothesis that the observed spectra are in 
approximate equilibrium. A distinct possibility is 
that, since normalized Boltzman rates for that cal- 
culation are 0(1), the character of the interaction 
is altered by non-resonant interactions. We work 
through this issue in the following four subsections. 
5.1.3.1. Representation Dependences 

The issue of representation (coordinate) depen- 
dent differences arises naturally in nonlinear prob- 
lems. Here it implies a question of whether a sta- 
tistically stationary state in one coordinate system 
similarly represents a statistically stationary state 
in another. 



Kinetic equations have been derived in Eule- 
rian, isopycnal and Lagrangian coordinate systems. 
Lvov et al. [2011] find that Eulerian, isopycnal and 
Lagrangian coordinate matrix elements in the hy- 
drostatic, non-rotating /^ <C o"^ ^ N'^ limit are 
equivalent on the resonant manifold. With ro- 
tation, comparisons are restricted to isopycnal co- 
ordinate formulation of Lvov and Tabak [2004] and 
the Lagrangian coordinate formulation of Olbers 
[1973]; McComas [1975]. Strikingly, the / / 0, hy- 
drostatic resonant expressions are also equivalent. 
Thus representation dependent differences are not 
apparent at the level of the RIA. 

Some physical sense can be made by noting that 
the resonant interaction approximation assumes, 
perforce, an expansion in terms of a non-advected 
wavefield, with linear dispersion relation (2). It is a 
description of the wavefield as a system of coupled 
oscillators with distinct resonance conditions. 

Non-resonant kinetic equations are limited to use 
of isopycnal coordinates. The isopycnal Hamiltoni- 
ans in Lvov and Tabak [2001] and Lvov and Tabak 
[2004] are explicitly expressed in terms of the the 
generalized coordinates and momenta. The gener- 
alized coordinates and momenta are implicit func- 
tions in the Lagrangian coordinate system formu- 
lations of Olbers [1976], McComas and Bretherton 
[1977] and Pomphrey et al. [1980] and the Eule- 
rian coordinate system formulation of [Voronovich, 
1979]. An implicit representation requires an ex- 
pansion in powers of small fluid parcel displace- 
ments in addition to an assumption of weak non- 
linearity. This small amplitude assumption repre- 
sents an unconstrained approximation whose do- 
main of validity vis-a-vis the assumption of weak 
nonlinearity is not well defined, [Miiller et al., 
1986]. This matters when the resonances are 
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broadened, Section 5.1.3.2. In contrast, the ex- 
plicit isopycnal coordinate representation of Lvov 
and Tabak [2004] requires only an assumption of 
weak nonlinearity. 

The next level of approximation to the ki- 
netic equation assumes wave amplitudes vary in 
time, a = a{T), and the energy conserving delta- 
functions in (37), d{ap — cjpi — o-pj), consequently 
need to be "broadened" to take near resonant in- 
teractions into account: When the resonant kinetic 
equation is derived, it is assumed that the ampli- 
tude of each plane wave is constant in time, or, in 
other words, that the lifetime of single plane wave 
is infinite. The resulting kinetic equation, never- 
theless, predicts that the wave amplitude changes. 
For small levels of nonlinearity this distinction is 
not significant, and resonant kinetic equation con- 
stitutes a self-consistent description. For larger 
values of nonlinearity this is no longer the case, and 
the wave lifetime is finite and amplitude changes 
need to be taken into account. A mathematical 
interpretation is that resonance broadening repre- 
sents a bandwidth Act connected to the temporal 
localization through a Fourier uncertainty principle 
, i.e., the product of wave lifetime r and bandwidth 
Ao" is given by rAo" = 27r. This effect is larger 
for stronger levels of the normalized Boltzman rate 
(48). The next section describes attempts to incor- 
porate resonance broadening effects into the kinetic 
equation. 

5.1.3.2. Renormalization of the Kinetic 
Equation 

Inclusion of resonance broadening requires a 
renormalization of the kinetic equation (37). De- 
Witt and Wright [1982, 1984] tackle this starting 
from the Lagrangian of Olbers [1976] and carry- 
ing forward their calculations in the context of 
the Direct Interaction Approximation (DIA) . Their 
nominally finite amplitude analysis still retains the 
small amplitude approximation inherent in the La- 
grangian formulation. Carnevale and Frederiksen 
[1983] begin from the Eulerian equations of mo- 
tion to derive the DIA closure equations for inter- 
nal waves in two dimensions and without rotation 
through the application of renormalization theo- 
ries originally formulated for quantum and classi- 
cal statistical field theory. A similar tack is pur- 
sued by Lvov et al. [1997] and a general repre- 
sentation for a 3 wave system is given in terms 
of a Hamiltonian representation. There the Wyld 
diagrammatic technique adapted for the statisti- 
cal description of a wave interaction system was 
used. This approach allows the Dyson- Wyld re- 
summation of the reducible infinite class of dia- 
grams, which presents certain parts of the nonlin- 
ear interactions as effectively being linear on aver- 
age. As a result, off-resonant interactions can effec- 



tively contribute to long-time statistical averages. 
This result is obtained by analytical resummation 
of the infinite diagrammatic series for the Green's 
function and double correlator. This result is given 
by 

^ = 4/ K,,pffpl2HP - PI - P2)d='pid3p2 
-4/ |l^pP/,pP/l2p-F(pi - P2 - p)d'pid^p2 
-4/ l^p':pj'/2pl-F(p2 - P - pOd^'pid^Ps , 



with -F(q) = (5q£(o-q). Here C is defined as 
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where Tki2 is the total broadening of each par- 
ticular resonance and Act represents the residual 
Ao" = o"d — Cd, — (Too, etc.. 
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The width of the resonance rfci2 in (60) is given 
by 

rfci2 = 7p+7pi +7p2- 
It means that the total resonance broadening is the 
sum of individual frequency broadening, and can 
be thus seen as the "triad interaction" frequency. 
The direct interaction approximation customarily 
sets 7p = 0. 

The single frequency renormalization is calcu- 
lated self-consistently from 

7p 

4 / I^Pi,P2l^("'Pi + '^P2)-^(P - Pi - P2)rfPlrfP2 

(np2 - npJJ'(pi - P2 - p)dpidp2 

-4 / l^p^piP("-Pi - "P2)-^(P2 - P - Pi)(ipi(ip2 

(61) 

Note that if the nonlinear frequency renormaliza- 
tion tends to zero, i.e., Tki2 — )■ 0, i2 reduces to the 
delta function: 

lim £(Acr) = 7r5(Acr). 
rfei2— >o 

In the limit of no broadening (59) reduces to (37). 
A self-consistent estimate of 7p requires the iter- 
ative solution of (59) and (61) over the entire field: 
the width of the resonance (61) depends on the 
lifetime of an individual wave [from (59)], which in 
turn depends on the width of the resonance (60). 
The result is a numerically intensive computation 
which has yet to be undertaken. Instead, Lvov 
et al. [2011] make the uncontrolled approximation 
that: 



7p 



ear 
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with e taken as 10^^ and larger. 

The patterns of gain and loss at e = 1 x 10^^ 
(Fig. 40) are quantitatively similar to previously 
published results including only resonant interac- 
tions, e.g., Olbers [1976]. Minimal variation is 
noted in evaluations using e = 1 x 10~^ and 
e = 1 X 10~^, implying that the e = 1 x 10~^ evalua- 
tions represent the resonant interaction limit. Nor- 
malized Boltzman rates are 0(1) at high wavenum- 
ber, and in particular within the high frequency 
part of the spectral domain which is nominally a 
stationary state in the ID approximation. [Lvov et 
al., 2011] find that normalized Boltzman rates are 
reduced by factors of 2-3 with e = 0.5 (See their 
Fig. 5). Here we present results for which it is dif- 
ficult to claim resonance broadening significantly 
decreases the normalized Boltzman rate (Fig. 40). 
The difference between results presented in Lvov et 
al. [2011] and in Fig. 40 is that Lvov et al. [2011] 
extend the GM76 spectrum to 2mc, zeroing the 
spectrum at higher vertical wavenumber. Here we 
roll off the GM76 spectrum as GM76xmc/m for 
m > rric, zeroing the spectrum only for m > Snic- 
We regard this as being a much more realistic char- 
acterization of the observations. Clearly, the nu- 
merical evaluations are sensitive to the treatment 
of the spectral boundaries. Experimentation sug- 
gests numerical evaluations of the broadened ki- 
netic equation are also sensitive to the functional 
representation of 7p other than (62). 

The iterative calculation needs to be carried out 
in order to address whether the patterns (43) asso- 
ciated with the Induced Diffusion mechanism (Fig. 
37) are more than simply fortuitous. 
5.1.3.3. Ray-tracing models 

The interpretation offered immediately above is 
one that attempts to define the breadth of the res- 
onance associated with a forced-damped system of 
weakly coupled oscillators, in which 'forcing' and 
'damping' can be construed as the scattering of 
energy into and out of (p,cr). An alternative de- 
scription is to view the wavefield as a collection 
of wave-packets, i.e., as particles. In this particle 
perspective, the time evolution of a wave packet 
(a 'test' wave) is proportional to gradients of the 
Doppler shift as one evaluates the evolution of the 
wavenumber, frequency and amplitude of the wave 
packet along its trajectory, in contrast to the cou- 
pled oscillator construct that assumes only an am- 
plitude modulation. These ray methods do not 
make a weak interaction assumption. They do, 
however, make scale separation and small ampli- 
tude assumptions. 

The ray equations [Witham, 1974] are: 



n 



dp 
lit 



for the time evolution of a wavevector following a 
wavepacket, and 



dr 
dt 



U + VryUJ 



(64) 



for the position of the wavepacket, which varies 
as the sum of advection plus group velocity. The 
background velocity u is assumed to be slowly 
varying in comparison to the wave phase and is 
treated as time-dependent stochastic realizations 
based upon the GM model in this wave-wave in- 
teraction problem. The intrinsic frequency 

uj = a — p u (65) 

is given by a linear dispersion relation (2). The am- 
plitude of the wavepacket is assumed to be small 
and obeys an action conservation principle: 



dn{p) 
dt 



+ Vr • (u + Cg)n(p) + Vp•7^n(p) = 0. 



(66) 

A rigorous derivation [Witham, 1974] is based upon 
the concept of phase (p • r — at) conservation for a 
single wave packet rather than energy (E?^]^o"jnj = 
0) and momentum (S^^^pjnj = 0) conservation of 
a triad system. The extent to which the dynamics 
inherent in the kinetic equation (59) are captured 
by the restricted description (63)-(64)-(66) has not 
been fully elucidated. 

The ray equations do not constitute a transport 
representation in (23), which attempts to describe 
the average evolution of a field, rather than the 
evolution of a single wavepacket. To this end, 
Flatte et al. [1985] report that test wave spec- 
tra exhibit a tendency to relax to one-dimensional 
vertical wavenumber spectra proportional to m~^ 
(GM76) and one-dimensional horizontal wavenum- 

1/2 

ber spectra proportional to /c^ times GM76, inde- 
pendent of whether the background was consistent 
with the GM model. This is the solution to the 
Kinetic Equation (37) defined by Pelinovsky and 
Raevsky [1977] and that which appears in the di- 
mensional analysis of Polzin [2004a]. The model 
diagnostics of Flatte et al. [1985] indicate action 
transport to higher frequency which is inconsistent 
with the induced diffusion model of action trans- 
port to higher vertical wavenumber without trans- 
port in horizontal wavenumber, i.e., transport to 
lower frequency. These tendencies are captured 
by a cascade representation of T^ found in Polzin 
[2004a]: 
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The cascade closure has been used in combination 
with PSI forcing of the inertial field [Polzin et al.. 
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Figure 40. Normalized Boltzman rates (48) for the GM76 spectrum calculated via (59) using the Lvov 
and Tabak [2004] matrix elements with S = 10"^ (left) and 5 = 5 x 10"^ (right), Lvov et al. [2011]. 
A value of 5 = 10""^ represents the resonant interaction limit. The spectra have been extended beyond 
m > rric by rolling off the spectrum as GMxrUc/m. 



2008] to obtain a link between power law depen- 
dencies in the vertical wavenumber and frequency 
domains. This closure is also used to discuss the 
near-boundary decay problem in Polzin [2004b] . It 
is equivalent to the ray-tracing based formulas em- 
ployed in diagnostic studies of dissipation parame- 
terizations [Gregg, 1989; Polzin et al., 1995; Gregg 
et al., 2003]. 

The most important point is the frequency do- 
main weighting of energy transports to small ver- 
tical scales in (67): 

.1/2 
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(68) 

Despite the small spectral density at high fre- 
quency, high frequencies make a significant con- 
tribution to the total spectral transports. 
5.1.3.4. Numerical Simulations 

Direct numerical simulations of the equations of 
motion are not limited by the dynamical assump- 
tions inherent in either the weakly nonlinear or 
eikonal representations. They are, however, sub- 
ject to computational restrictions. 

Winters and D'Asaro [1997] present spin-down 
simulations based upon the G1VI76 spectrum with 
varying amplitude. They regard their results as be- 
ing consistent with the magnitude (e = 7 x 10"^'' 
W/kg) and scaling of Gregg [1989], e ~ E'^N'^. 
Polzin et al. [1995] and Polzin [2004a] argue that 
theoretical estimates based upon the eikonal model 
of [Henyey et al., 1986] and the dynamic balance 
of JVIcComas and IVIiiller [1981b] are indistinguish- 
able from the observational results. Thus, both 
observations, Henyey et al. [1986], McComas and 
Miiher [1981b] and Winters and D'Asaro [1997] 
are consistent. The domain considered by Win- 



ters and D'Asaro [1997] consists of a rectangular 
box 80 X 10 X 1 km on a side with resolved wave- 
lengths of 1 km in the horizontal and 50 m in the 
vertical. Note that this domain does not include 
regions in Fig. 38 exhibiting large values of the 
normalized Boltzman rate. Interactions in the re- 
solved domain may be dominated by PSI transfers 
as discussed in IVIcComas and Miiller [1981b]. 

Forced non-rotating simulations are presented 
in Furue [2003]. The computational domain is a 
box of horizontal size 100 x 100 x 128 m height. 
The forcing is isotropic in wavenumber and peaks 
at a horizontal wavelength of 25 m. The forc- 
ing is controlled so that amplitudes are consis- 
tent with GM76 and the resulting dissipation is a 
significant fraction of that associated with GM76. 
While we are uncertain about how to relate these 
forced simulations to the wave-wave interaction 
problem, we believe finding significant dissipation 
at high frequencies is credible [see (67) and Section 
5.4.2]. Also notable are findings that wavenumber- 
local interactions dominate the transfers and that 
energy transport to higher horizontal wavenum- 
bers (Section 5.1.1.e) is a robust feature of the 
simulations whereas energy transport in vertical 
wavenumber is strongly dependent upon spectral 
shape. 

Quantitative estimates of spectral transports in 
high resolution and realistically forced direct nu- 
merical simulations is a high priority, especially in 
conjunction with efforts to understand the limita- 
tions of a broadened kinetic equation (59). 
5.1.4. Does the Continuum Exist? 

Eulerian frequency spectra at high vertical 
wavenumber are contaminated by vertical Doppler 
shifting: near-inertial frequency energy is Doppler 
shifted to higher frequency at approximately the 
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Figure 41. Fig. 4b from Pinkel [2008] witli a dis- 
placement variance overlay and associated velocity spec- 
trum, assuming the Gaussian distribution given by Pinkel 
[2008] with 2 hour lag correlation and 1.25 m^ displace- 
ment variance. The associated velocity spectrum (upper 
curve within the overlay) lies above the observed velocity 
spectrum. A buoyancy frequency oi N = 6 cph has been 
assumed. 

same vertical wavelength. High frequency, small 
horizontal scale internal waves are also easily ad- 
vected past a similar Eulerian observer by the hor- 
izontal velocity of near-inertial waves. Such differ- 
ences motivate our concern with different coordi- 
nate systems in Section 5.1.3.1. 

Pinkel [2008] posits that the apparently continu- 
ous wavenumber-frequency spectrum at high verti- 
cal wavenumber results from the advective smear- 
ing of variance contained in a few discrete spectral 
lines. A "simple" model is introduced to quan- 
tify the effects of lateral advection, and both ran- 
dom and deterministic vertical advection on these 
spectral lines. Two data sets, including one from 
the Arctic (the Surface Heat Budget of the Arctic 
[SHEBA], Fig. 41) were used to demonstrate the 
model. 

Model implementation includes arbitrary choices 
of functional forms for lag coherence estimates and 
choosing model parameters to maximize the high 
frequency shear variance which would be inter- 
preted as resulting from Doppler smearing. Fur- 
ther observational constraints, such as requiring 
the vertical displacement spectrum to be consis- 
tent with the observations, could have been im- 
plemented. As presented, the model's vertical dis- 
placement spectrum appears about a factor of two 
larger than expected (Fig. 41). 

The SHEBA spectrum is quite different from the 
AIWEX spectrum. The SHEBA frequency domain 
power laws are much steeper. This is possibly re- 
lated to time dependence in forcing or nonlinearity 
(the SHEBA spectrum was obtained in August- 
September 1998, AIWEX was conducted in March- 
May 1985 and Pinkel [2005] notes significant tem- 
poral variability in the stratification with the re- 
cent buildup of a seasonal thermocline prior to the 



SHEBA observations) and/or to spatial variability 
(the reported SHEBA observations were obtained 
atop the continental shelf whereas the AIWEX 
data represent deep ocean conditions). The AI- 
WEX data Levine [1990] appear to argue strongly 
for a high frequency, high wavenumber continuum. 

While there assuredly is some contamination 
of the high vertical wavenumber continuum by 
Doppler shifting, the goodness of fit in Pinkel 
[2008] can not be interpreted as evidence for the 
nonexistence of the continuum. If robust, the 
goodness of fit simply implies that the velocity ob- 
servations are consistent with the proposed model. 
Further consistency tests based upon both den- 
sity and velocity fields could demonstrate other- 
wise. Note that the 2-D PATCHEX and SWAPP 
spectral parameters (Section 3.3.5) were extracted 
from Ep rather than E^ spectra in an attempt to 
limit the effects of Doppler shifting and instrument 
response functions. 
5.1.5. A Summary and Future Challenges 

The characterization of nonlinear interactions 
is a complicated process and therefore reductions 
to simplified forms for various limiting cases are 
highly desirable. The reduction of the kinetic equa- 
tion in the ID limit to a Fokker-Plank equation in 
vertical wavenumber is fundamental in this regard: 
the stationary states of this equation collapse much 
of the observed variability in the spectral parame- 
ters (a, b) (Fig. 37). There are, however, significant 
questions regarding the validity of the underlying 
assumption of resonant interactions, especially in 
the ID limit at small vertical scales, that prompt 
one to wonder if the pattern match is merely for- 
tuitous. At these small vertical scales, diagnostic 
studies that compare simplified representations of 
the downscale transport with estimates of turbu- 
lent dissipation often settle upon the ray-tracing 
model (67) as providing the theoretical framework 
for the finescale internal wavefield. However, the 
simplified representations derived within the RIA 
provide predictions that are similar both in mag- 
nitude and in scaling [Polzin, 2004b]. Thus one 
might be prompted to suggest that such misgiv- 
ings are simply due to offenses against one's theo- 
retical sensibilities. 

There is, however, a very substantive issue that 
cries out for explanation. Associated with the ID 
mechanism and b > are energy transports to 
higher vertical wavenumber and lower frequency. 
A stationary state under the ID mechanism and 
b > thus requires either: 

• An external energy source at high frequencies 
External forcing is discussed in Section 5.2 and 5.3 

• An alteration of spectral transports associ- 
ated with finite amplitude effects that acts to re- 
duce transports through the frequency domain. 
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Section 5.1.3 was dedicated to the issue of reso- 
nance broadening. 

• A reassessment of resonant transfers. We find 
that residuals in the resonant calculation for GM76 
give rise to normalized Boltzman rates ep that are 
0(1). We argue that these residuals arise from 
inherently local interactions and hypothesize that 
they are associated with spectral transports in the 
horizontal wavenumber domain. It is possible that 
such residual interactions determine the transport 
of GM76. 

• An effective source associated with nonhydro- 
static effects. Only the numerical evaluations of 
Olbers [1976] have included nonhydrostatic effects. 
Those results are silent on this issue. 

5.2. Wave-Mean Interactions 
5.2.1. Spectral Patterns 

Much of the analysis in Section 5.1 was predi- 
cated upon the assumption of an inertial subrange; 
that there exists a range of wavenumbers and fre- 
quencies removed from the influence of forcing and 
dissipation, and far from the characteristic spatial 
and temporal scales of the system, that the spec- 
trum can be characterized as being scale invari- 
ant. We have noted the difficulties for this sce- 
nario with respect to the importance of extreme 
scale separated interactions in Section 5.1.1. Here 
we note that wave-mean interactions pose a spe- 
cial problem for the scale invariant solution in that 
wave-mean interactions can be cast as a source rep- 
resentation So{p) acting throughout the spectral 
and spatial domains rather than being located at 
the boundaries. 

A large body of work focusses upon the simpli- 
fied cases of waves interacting with parallel shear 
flows [Molemaker et al., 2005; Ruddick, 1980] or 
vortices [Kunze, 1985], i.e., symmetric background 
flows. This conceptual paradigm may be made 
much more general and useful by allowing the 
background flow to have spatial gradients in both 
horizontal coordinates. While the assumption of 
a spatially symmetric background flow makes for 
a more tractable theoretical problem, it also in- 
troduces a hidden constraint in which the verti- 
cal flux of streamwise pseudo-momentum is non- 
divergent except at a critical layer. In the gen- 
eral problem of internal waves interacting with a 
mesoscale eddy fleld having 3-dimensional struc- 
ture, a wave packet can freely exchange wave 
psuedo-momentum for eddy potential vorticity, 
Biihler and Mclntyre [2005] . This is restricted from 
happening in a symmetric background flow. 

The essence of the interaction in three dimen- 
sions is the filamentation of an internal wave by 
the horizontal rate of strain and vertical propaga- 
tion into a critical layer. This can be directly seen 



from a linearized internal wave energy equation 
for wave propagation in a quasigeostrophic back- 
ground [Polzin, 2010]: 

d 



(q-^+u- ^h){Ek + Ep)+V- tt"u" = 



4-2^"v"P^ 



[v"w" + -^V'u"]v, (69) 



with kinetic [Ek = (u"^ + v"^ + w"^)/2] and po- 
tential [Ep = {N~'^b"'^)/2] energies. Temporal 
variability and advection of internal wave energy 
by the geostrophic velocity field are balanced by 
wave propagation and energy exchanges between 
the quasigeostrophic and internal wave fields. Non- 
linearity and dissipation are assumed to be higher 
order effects. Spatial gradients of the vertical ve- 
locity w do not appear as Wz is small [order Rossby 
number squared] in the quasigeostrophic approxi- 
mation. The thermal wind relation has been in- 
voked to cast the vertical Reynolds stress and hor- 
izontal buoyancy flux as the rate of work by an ef- 
fective vertical stress acting on the vertical gradient 
of horizontal momentum. These two terms in the 
effective stress will cancel each other in the limit 
that oj — > /, Ruddick and Joyce [1979]. Finally, 
Ss = Vx + Uy represents the shear component and 
Sn = Ux — Vy the normal component of the rate of 
strain tensor. In the limit of an extreme scale sep- 
aration between wave and background flow, hor- 
izontal group velocities are small and an internal 
wave haves much like a passive tracer, [Biihler and 
Mclntyre, 2005; Polzin, 2008a] as the rate of strain 
cascades the wave to higher horizontal wavenum- 
ber. As it does so, the horizontal components of 
the internal wave stress result in work against the 
mean flow gradients and gain of wave energy. In 
the situation in which a wave is propagating verti- 
cally into a critical layer, the vertical stresses are 
such that the wave looses energy to the eddy field. 
The coupling of mesoscale eddies and internal 
waves and construction of a So{p) representation 
was explicitly addressed in Miiller [1976]. Miiller 
[1976] 's analysis was based upon a perturbation 
expansion of the internal wave radiation balance 
equation (23). Internal waves were assumed to be 
of small amplitude and have small spatial scales 
relative to a geostrophically balanced background 
that evolves over a much longer time scale. Spatial 
gradients in the background were assumed to be 
sufficiently weak that wave-mean interactions af- 
fect wave propagation only through an advective 
Doppler shift. In this limit of small amplitude 
waves interacting with a quasi-qeostrophic eddy 
field, exchanges of energy, momentum and vortic- 
ity are reversible unless an external force is invoked 
which acts in phase on the wave-mean induced per- 
turbations. In Miiller [1976], the external force is 
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identified as the tendency for nonlinear interactions 
to relax the wave-mean induced perturbations to 
an isotropic state. It is this relaxation that cre- 
ates a permanent exchange of pseudo-momentum 
for potential vorticity. In order to close the prob- 
lem, Miiller [1976] invoked a relaxation time ap- 
proximation and further assumed the relaxation 
time scale to be constant. Estimates of the eddy- 
internal wave coupling strength were not supported 
by observations [Ruddick and Joyce, 1979; Brown 
and Owens, 1981] and this line of investigation was 
dropped. With regards to the strawman in Fig. 34, 
the consequence is that wave-mean interactions are 
not available as a source to supply energy depleted 
by PSI transfers. 

Miiller [1976] was written from the perspective 
that the thermocline is characterized by a diapy- 
cnal diffusivity of Kp = 1 x 10~^ m^ s~^, which 
implies a time scale of 5-10 days for nonlinear in- 
teractions to drain energy out of the background 
internal wavefield. Several decades of research has 
since demonstrated that the background diffusiv- 
ity is Kp = 5 X 10^^ m^ s^^, with corresponding 
time scale of 50-100 days. The difference is crucial, 
as Miiller assumed that the eddy-internal wave in- 
teraction was local in the spatial domain. With 
0(50-100) day relaxation time scales, wave prop- 
agation effects become important: larger scale in- 
ternal waves can propagate through an eddy-wave 
interaction event with minimal permanent ex- 
change of pseudo-momentum and vorticity. For- 
tunately, Miiller [1976] tells us how to account for 
such propagation effects. Accounting for propaga- 
tion effects and use of a more realistic prescription 
of the nonlinear time scale returns estimates of the 
coupling that are in much better agreement with 
the observations, [Polzin, 2008c]. 

The radiation balance equation source function: 

aSo{v) = S^ + S: (70) 

can be characterized with the following relaxation 
time scale formulas: 
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with relaxation time scale tr and propagation time 
scale (wave group velocity divided by eddy length 
scale) Tp. Evaluations of the relaxation time-scale 
formulas are presented in Fig. (42). We use 
two spectra intended to mimic those at Site-D 
(GM76) and the Sargasso Sea [(18) and (19) with 



(s = 2,t = 2.3,i* = 15,r = 0.8, q = 0.5, iV = 
2.6cph,^ = 0.7^^^^^)]. Mean flow gradient statis- 
tics are taken from the LDE data [Polzin, 2010]: 
{H, L, ul + vl, ul + ul + vl + ij2)=(700 m, 100 
km, 4.1 X 10^^ s^^, 5.8 x 10^-*^^ s^^). Preliminary 
estimates from Site-D array data indicate that this 
is an appropriate characterization there. The re- 
laxation time scale is specified using (67). 

The coupling is a complicated function of the 
spectral parameters: spectral amplitude, high 
wavenumber/frequency power law specifications 
and low wavenumber roll-off; and external param- 
eters, i.e., mesoscale velocity gradient variances 
and characteristic spatial scales (Fig.s 45 and 44). 
Without having executed a complete parameter 
regime survey, we rationalize the results as follows. 

In the continuum frequency band, differences be- 
tween the vertical source functions for GM76 and 
the Sargasso Sea result from the GM76 3-d action 
spectrum having minimal gradients of action in the 
vertical wavenumber domain, i.e., 6 = 0. For b ^ 0, 
the zero-contour demarcating source/sink tenden- 
cies is dictated by m*: larger vertical wavelengths 
are sources of energy for the internal wavefield. 
The magnitude of the positive region relative to 
the high wavenumber negative region is dictated by 
the ratio of relaxation to propagation time scales, 
Tr'-Tp. The possibility of a positive feedback mech- 
anism dictating m* is obvious. 

In the near-inertial band, vertical transfers are 
sinks and horizontal transfers are sources, consis- 
tent with a vertical critical layer and horizontal fil- 
amentation scenario. The net transfer obviously is 
a function of the horizontal vs vertical gradient 
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variances of the mesoscale eddy field, i.e., it's as- 
pect ratio. The Sargasso Sea estimates imply an 
eddy aspect ratio slightly less than f /N and this 
may be typical [Arbic and Flierl, 2004]. We have 
also included non-hydrostatic conditions in these 
calculations. There is a signature of a significant 



source at the buoyancy frequency in the vertical 
source function that is largest for 6 7^ 0. This 
could potentially resolve issues related to ID trans- 
fers requiring a high frequency source in order to 
maintain a transport of action to high vertical 
wavenumber. Section S.l.l.d. Some skepticism 
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Figure 42. The spectral source functions defined in (71) and (72). Left-hand panels are GM76 eval- 
uations. Source functions in the right-hand panels use a 'Sargasso Sea' spectrum. Panels A, B, and C 
represent the horizontal, vertical and combined source functions. The thick line separates positive source 
values from negative. Negative values of the source function are, effectively, energy sinks for the internal 
wavefield. 
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Figure 43. The spectral source functions appearing in Fig. 42 but expressed as a normalized trans- 
fer rate 2nSo{o',m)/aE{a,m). This normalization permits a direct comparison with Fig.s 38 and 40. 
Contour intervals here are 0.01. 
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is required: such dynamics will compete with turn- 
ing point effects and questions of how to specify H 
in Tp arise. Apart from the buoyancy frequency, 
note that these wave-mean transfers are far too 
small to offset resonant wave- wave interactions: 
Compare Fig. 38 with 43. 
5.2.2. Budgets 

Energy exchange between the mean and wave fields 
results from a stress-strain relation in the hori- 
zontal. Direct evaluation of the energy exchange % 
using the PolyMode LDE current meter array ■% 
data [Polzin, 2010] returns a transfer rate some- -; 
what smaller than previously obtained [Brown and 
Owens, 1981], 4 x 10" 1° W/kg at 825 m depth. 
If this value characterizes the thermocline region 
having a characteristic depth H = 1000 m, a net 
source strength of 0.4 mW m~^ can be anticipated. 
Table 1. This estimate may be biased low, how- 
ever, as both vertical and horizontal velocity gra- 
dients increase towards the surface and the rate 
at which energy is transfered is proportional to 
the velocity gradient variances. The LDE was lo- 
cated at (31°N, 70° W), midway between Lotus and 
Fasinex. The estimates of energy transfer through 
eddy-internal wave coupling at 825 m depth do not 
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appear to explain the regional seasonal cycle of in- 
ternal wave energy. There is, however, a seasonal 
cycle in upper ocean fronts associated with buoy- 
ancy forcing, e.g., Weller [1991]. This upper ocean 
frontal regime is characterized by an 0(1) Froude 50' 
(vertical gradient of horizontal velocity divided by 
buoyancy frequency) and Rossby (relative vorticity 
divided by /) number parameter regime [Pollard 
and Regier, 1992] which may support more effi- 
cient energy transfers than the quasi-geostrophic 
regime at the thermocline base. 



40 



3 



5.3. Forcing 

The focus of the previous subsections was the gen- ^ 
esis of spectral structure. In Sections 5.3, 5.4 and 
5.5 we focus more upon (i) whether the variabil- 
ity in spectral parameters represents the character 
of the forcing at the endpoints of the cascade / 
spectral boundaries, and (ii) whether the forcing 
is consistent with the direction and magnitude of 2(f 
energy transfers. 

5.3.1. Near-inertial wave generation by a 
variable wind stress 0^ 

Work done by the wind in generating inertial fre- 
quency mixed-layer motions has been the subject Figure 46. Map of near-inertial wind input into the 
of several recent publications. Nearly- global maps mi^ed layer in units of W/m^ [Alford, 2003]. 
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Figure 48. Power lost from the internal wavefield to the surface wavefield in units of 10 '^ W/m^ for 
summer months [Watson, 1994]. 
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Figure 49. Power lost from the internal wavefield to the surface wavefield in units of 10 '^ W/m^ for 
winter months [Watson, 1994]. 



of this vifind work using the Pollard and Millard [1970] 
damped-slab model are presented in Alford [2003]; 
Watanabe and Hibiya [2002], Fig.s 46 and 47. The 
slab model is highly idealized. Using observations of 
the upper ocean response to atmospheric forcing and 
a more sophisticated mixed layer model, Plueddemann 
and Farrar [2006] describe how the slab model tends 
to overestimate inertial frequency wind work by a fac- 
tor of two. Their estimates for the Lotus, Fasinex and 
Subduction field programs are given in Table 1. Similar 
estimates for Site-D are presented in Silverthornc and 
Toole [2009]. 

Plueddemann and Farrar [2006] also note that other 
modeling investigations (Crawford and Large [1996]; 
Skyllingstad et al. [2000]) imply that only a fraction 
of this work (50%) may be available for near-inertial 



internal wave generation as the other part is lost to 
viscous dissipation in the mixed layer for the resonant 
forcing conditions that dominate the revised transfer es- 
timates. 'Best' estimates of source strength in Table 1 
incorporates this dissipation. 

The aforementioned work concerns only getting inertial 
energy into the mixed layer. Radiation of this energy 
into the stratified interior can be represented by invok- 
ing a Fourier decomposition of the slab model and map- 
ping onto the internal wavefield using a linear model 
that includes the variation of the Coriolis parameter 
with latitude [D'Asaro et al., 1995]. There are, however, 
unresolved issues: the linear model under estimates the 
initial rate at which energy leaves the mixed layer by 
20 — 50% and can not explain the evolution of the small 
vertical wavelength 
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near-inertial field. Nonlinear interactions are held to 
be a primary candidate to explain the later [D'Asaro, 
1995]. For further consideration of near-inertial motions 
see Section 5.5. 

5.3.2. Surface-internal wave coupling 
Internal waves can also exchange energy with the sur- 
face wavefield, thereby implying a parametric depen- 
dence of source strength upon wind stress and possible 
seasonality. In the resonant interaction approximation, 
a pair of surface waves with nearly equal wavenumber 
and frequency can interact with an internal wave mode 
through a difference interaction, Watson et al. [1976]; 
Olbers and Herterich [1978]. Estimates of energy trans- 
fer are sensitive to the structure of the upper ocean 
buoyancy profile. Estimates [Watson, 1994] using cli- 
matological buoyancy profiles and seasonally averaged 
wind estimates return transfer rates of several tenths 
of a mW m^^ from the internal wavefield to the sur- 
face wavefield during summer months (Fig. 48) and at 
low latitudes during winter months (Fig 49). The high 
latitude transfer rates during winter are smaller. In 
contrast to these wind wave results, rapid transfers of 
energy from a narrow band ocean swell to the in- 
ternal wavefield are possible [Watson, 1990] and remain 
to be quantified. 

The coupling of surface and internal gravity waves 
thus does not appear to be the effective source at high 
frequencies and low vertical wavenumber as depicted in 
the strawman Fig. 34 and required as a boundary con- 
dition on the energy transports in the dynamic balance. 
Section 5.1.1.d. 

5.3.3. Internal wave coupling to mixed- 
layer turbulence 

The coupling of surface and internal gravity waves is 
not the only possible source of high frequency internal 
wave energy. It has long been conjectured that mixed- 
layer turbulence can be coupled to the internal wave- 
field. Bell [1978]. Observations presented in Wijesekera 
and Dillon [1991] reveal a narrow band internal wave 
response in the upper equatorial ocean at frequencies 
near the upper ocean stratification rate. Skyllingstad 
and Denbo [1994] point towards shear instability as a 
generation mechanism. It is unclear how this mecha- 
nism translates to a rotating system. 

Polton et al. [2008] use large-eddy simulations of 
mixed-layer turbulence with a parameterization of 
Stokes drift forcing by a surface gravity wavefield. This 
forcing sets up Langmuir circulations and a inertial cur- 
rent oscillation is superimposed that advects the Lang- 
muir cells relative to the stratified interior. Vertical 
circulations of the Langmuir cells excite high frequency 
internal waves having a phase velocity that matches the 
depth averaged mixed layer velocity. 

Polton et al. [2008] estimate that energy and momen- 
tum are radiated at rates comparable to those associ- 
ated with near-inertial wave radiation. These waves, 
however, are likely trapped within the highly stratified 



transition layer with either implications for mixing or 
for acting as a source of high frequency energy that can 
be transported to lower frequency via wave-wave inter- 
actions. Analysis of upper ocean data sets for a high 
frequency response on the time scale of the wind forcing 
could be quite enlightening. 
5.3.4. Tides 

5.3.4.1. Spatial Patterns 

Numerical models have, in several instances, been 
used to both estimate the energy conversion from 
barotropic to baroclinic tides and keep track of the 
nonlinear evolution of the internal tides (e.g., Simmons 
[2008]). While such models may not have the spatial 
resolution required to get the details of the nonlinear 
transfer process right, or pick up the generation of baro- 
clinic tides from small scale topographic roughness, they 
should give some indication of where the PSI decay can 
expected to be important. The eastern Atlantic is one 
such region, Fig. 50, providing further circumstantial 
evidence that the Natre spectrum is forced primarily 
through PSI decay. 

5.3.4.2. Budgets 

The issue of internal tide generation, propagation and 
decay has received a great deal of attention this past 
decade, in part due to several major field programs 
(the Brazil Basin Tracer Release Experiment Ledwell et 
al. [2000] and the Hawaiian Ocean Mixing Experiment 
Rudnick et al. [2003]). Near global maps of the rate of 
conversion of barotropic tidal energy to baroclinic (in- 
ternal) tidal energy have been provided by several Eg- 
bert and Ray [2003]; Simmons et al. [2004]; Simmons 
[2008]; Niwa and Hibiya [2001] groups. As our focus 
is on the background internal wavefield, none of the 
data sets described here is from a region of significant 
barotropic to baroclinic conversion. 

The relevant phenomena appears to be the far-field 
decay of a low-mode tide via resonant interactions. This 
particular phenomena is believed to be most impor- 
tant when the frequency of the low-mode tide is slightly 
higher than the Coriolis frequency, e.g., Furuichi et al. 
[2005]. For a ~ M2, this corresponds to latitudes equa- 
torward of 28.9°. Many of the data sets examined here 
are in the vicinity of this critical latitude, particularly 
in the western North Atlantic. There, however, esti- 
mates suggest a low mode internal tide with horizon- 
tal energy flux of approximately 300 W m~^ [Hendry, 
1977]^ and analyses of relevant data from the western 
North Atlantic have not emphasized an interpretation 
of semidiurnal tidal PSI decay. An order of magni- 
tude estimate of source strength as the observed flux 
distributed over the distance from the assumed source 
(the Blake Escarpment some 700 km distant) is smaller 
( 300 W m^^ / 700 km) than the nominal background 
dissipation rate, 0(1 mW m~^). Section 5.4. Moreover, 
the internal tide likely cannot explain the observed sea- 
sonal cycle in high frequency internal wave energy if the 
source has been correctly identified. Baines [1982] finds 
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Figure 50. Maps of baroclinic tidal and near-inertial energy. From Simmons [2008]. 



that the details of the stratification profile are impor- 
tant in estimating the rate of conversion of barotropic 
tidal to baroclinic energy, but the Blake Escarpment is 
a relatively deep feature and extends only to 1000 m 
water depth, well below seasonal changes in stratifica- 
tion. 

Decay of a low mode-tide via resonant interactions 
poleward of the critical latitude would be through the 
induced diffusion mechanism. Olbers and Pomphrey 
[1981] argue that decay of a line spectrum with a typ- 
ical energy density of 100 J m~^ would result in loss 
from the tides at rates between 0.005 to 0.05 mW m~^. 
Such values will not contribute significantly to the in- 
ternal wave energy balance, but note; the calculation 
is based upon GM76, and with 6 = the estimate is 
biased. 

Barotropic tidal velocities in the eastern North At- 
lantic are typically larger than in the west, and, con- 
sistently, available estimates of the low mode internal 
tide (Rossby [1988]; Siedler and Paul [1991]; Alford and 
Zhao [2007]) are typically larger. However, such esti- 
mates exhibit significant spatial and temporal variabil- 
ity. Preliminary estimates of internal tide parameters 
from the Subduction array's center mooring assign rms 
buoyancy scaled velocities Uo, Vo of (0.02,0.02) m s^^ 
and 4 m displacement to mode-1, and implying an en- 
ergy fiux of approximately 1 kW m^^. This is a surface 



mooring and contamination of energy flux estimates by 
vertical motion in response to ambient currents is a def- 
inite concern. Alford and Zhao [2007] present energy 
fiux estimates of roughly 750 W/m for this mooring. 
This energy fiux is directed to the northeast, seemingly 
at odds with both near-field and far-field sources (a mi- 
nor seamount that rises to depths of 3000 m immedi- 
ately to the east, the Great Meteor seamount to the 
north and the European continental shelf to the north- 
east). Energy flux vectors in Simmons [2008] are to the 
west. 

An order of magnitude estimate of the flux divergence 
(Table 1) can be obtained by dividing 1 kW m~^ by the 
length scale characterizing the spatial variability in the 
internal tide amplitude, roughly 1,000 km. 

5.4. Dissipation 

Our underlying hypothesis is that the observed vari- 
ability can be understood as spatially local stationary 
states of the radiation balance equation (23). In this 
context, the role of dissipation is possibly two fold. 
First, in a local stationary balance, forcing balances dis- 
sipation. Second, dissipation may play a role in actively 
shaping high wavenumber spectra. We take up these is- 
sues in turn below. 
5.4.1. Budgets 
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TABLE 1. Estimates of source strengths. Best estimates of the internal wave source 
associated with wind work are a fraction (50 %) of the wind w^ork. The internal wave 
source associated with the decay of a low-mode tide is obtained by dividing the energy flux 
estimate by a length scale. Those length scales are ill constrained. 



program 


latitude 


longitude 


mechanism 


start 


stop 


source strength 


Lotus 


34° N 


70°W 


wind work 


May 1982 


Oct. 1982 


1 mW m'-'x 50% 


LDE 


31°N 


70°W 


eddy-wave coupling 


May 1978 


July 1979 


> 0.4 mW m-2 


Mode 


28°N 


70°W 


tide 


March 1973 


July 1973 


300 W m-i / 700 km 


Fasinex 


27°N 


70°W 


wind work 


Jan. 1986 


May 1986 


0.4mWm-2x 50% 


Subduction 


25.5°N 


29°W 


wind work 


April 1992 


Oct. 1992 


0.05 mWm-^x 50% 


Subduction 


25.5°N 


29°W 


tides 


April 1992 


Oct. 1992 


1 kW m"i / 1000 km 


Site-D 


39°N 


69°W 


wind work 


November 2001 


April 2006 


1.0 mW m-^x 50% 



Internal wave energy can be dissipated in the interior 
through wave breaking and at the boundaries as work 
done against viscous stresses. 

Depth integrated dissipation rates in the open ocean 
win typically exceed dissipation in the bottom bound- 
ary layer, D'Asaro [1982]; Fu [1981]. The contribution 
of the bottom boundary layer is mitigated by buoyancy 
scaling of internal wave velocity and low stratification 
in the abyss. Away from regions of topographic genera- 
tion, dissipation associated with interior wave breaking 
is dominated by thermocline and upper ocean regions 
having large stratification rates. One atypical exam- 
ple in which boundary layer dissipation can dominate 
is under ice cover in the Canadian Basin of the Arc- 
tic Ocean. There, the internal wavefield is anomalously 
weak and buoyancy scaling implies internal wave veloc- 
ities are largest near the upper boundary. A second 
atypical example is in regions of large geostrophic cur- 
rents in which wave dissipation can be modeled through 
a linearization of the bottom stress [Polzin, 2010]. 

Interior wave breaking is usually conceived of as a 
high wavenumber sink associated with instability mech- 
anisms, e.g., Polzin [1996]. That sink represents a pro- 
duction term 7^ in a turbulent mechanical energy bud- 
get: 



aS,=V = B- 



(73) 



being balanced by a buoyancy flux B and dissipation e 
and a stochastic production/dissipation balance. The 
correlations inherent in both production and buoyancy 
flux terms are difficult to measure directly. Rather, spe- 
cialized instrumentation and platforms are used to es- 
timate the rate of dissipation e. Full-depth estimates of 
such microstructure data are limited. In its stead, use 
has been made of flnestructure parameterizations based 
upon (67). These assign a value of 

8 X 10^1° W/kg 



V = 
and 
Kp = B/N^ = 5 X 10"*^ 



m^s-i 



to GM76 at 30 degrees latitude and a buoyancy fre- 
quency of 3 cph. 

With either fine- or microstructure method, a major 
caveat is that one needs to understand how (not just 
whether) the sampling of the wavefield is biased. For 
example, while one might hope to make sense of obser- 



vations in NATRE as being almost exclusively forced 
by a PSI decay [Polzin et al., 2008], analysis of the di- 
apycnal dispersion of an anthropogenic tracer suggests 
a long term average of a factor of two larger. The differ- 
ence may actually be related to the spatial variability 
of the tide. Experience [Polzin, 2009] suggests that is- 
sues of biasing in regions of rough topography are much 
more problematic than the factor of two in NATRE. 
5.4.2. Spectral Patterns: models 

In order to reconcile the 0(1) high wavenumber nor- 
malized Boltzman rates with the claim that GM76 rep- 
resents a stationary state, McComas and Miiller [1981b] 
invoke a spectral representation that assumes dissipa- 
tion events have space-time scales much smaller than 
the scales of the internal waves and that the turbulence 
can be parameterized by an eddy viscosity v, resulting 
in 



^S,(p) = 



dEjp) 
dt 



'E{p). 



(74) 



Similar representations are examined in greater detail 
in Natarov and Miiller [2005]. 

Some generalization of this model may be required if 
it is necessary to invoke dissipation in order to obtain a 
stationary state. Dissipation events will typically have 
vertical scales smaller than the nominal 10 meter cut- 
off, gauged either by correlation length scales of turbu- 
lent dissipation [Gregg et al., 1993] or the exclusive con- 
tribution of sub- 10 meter wavelength shear to the occur- 
rence of supercritical Richardson number events [Polzin, 
1996] (with supercritical Richardson number being in- 
dicative of shear instability), so the spatial scale sepa- 
ration between waves and turbulence could be justified. 
But the application of phenomenological mixing models 
indicates that dissipation events have lifetimes in excess 
of a buoyancy period [Polzin, 1996] so that a tempo- 
ral scale separation is problematic. In general, Polzin 
[1996] finds such phenomenological mixing models to 
be useful diagnostic tools, but difficult to develop into 
prognostic models because of this finite lifetime. 

A key point in a generalized model is that the 
flnestructure parameterization (67) weights the contri- 
bution of each frequency to dissipation as a product 
of vertical gradient variance and aspect ratio and thus 
the contribution of high frequency waves. The dissipa- 
tion function (74) weights each frequency only as the 
vertical gradient variance and so emphasizes the near- 
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inertial that dominate the vertical shear. Phenomcno- 
logical mixing models [Polzin, 1996] have been shovi^n 
to be in agreement with (67). Observations presented 
in Alford and Pinkel [2000] also appear to support the 
role of high frequency variability. 

A qualitative distinction between (74) and (67) may 
be drawn by noting that (74) is likely to be appropriate 
for a system of coupled oscillators. It may not make 
sense for a description of the internal wavefield as a 
system of wavepackets in which dissipation events are 
associated with the transport of individual wave packets 
to breaking scales, as invoked in the ray tracing mod- 
els, e.g., Henyey et al. [1986]. The ray-tracing model 
does not need a dissipation function because the wave 
packets are followed as particles to their annihilation. 

We regard the impact of dissipation on the internal 
wave spectral domain as an open question awaiting rig- 
orous analysis. 

5.4.3. Spectral Patterns: The mid-frequency 
dip and bottom boundary conditions 

Many of the frequency spectra are poorly character- 
ized by a simple power law fit near semidiurnal frequen- 
cies (e.g.. Fig. 14). Apart from the tidal peak, observed 
spectral levels tend to be smaller than the power law 
characterization. Levine [2002] forwards the hypothesis 
that this mid-frequency dip is related to nonlinearity in 
how energy is transfered to higher and lower frequency 
out of the semidiurnal peak. Here we forward an ad- 
ditional hypothesis regarding such departures. With 
a buoyancy profile that decreases towards the bottom 
boundary, only the low frequency portion of the wave- 
field reaches the bottom, where wave energy can be 
dissipated as work done against viscous stresses and 
vertical scale transformations associated with reflection 
or scattering from a nonuniform bathymetry can act. 
Calculations have already been made to assess the rate 
at which energy is redistributed by such scale transfor- 
mations for the GM spectrum, e.g., [Miiller and Xu, 
1992; Eriksen, 1985b]. With an eye toward explaining 
spatial variability of the background spectrum, one can 
also ask the question of how the equilibrated spectrum 
could be shaped by the cumulative effects of scale trans- 
formations acting through the bottom boundary condi- 
tion and nonlinearity acting through radiation balance 
equation's source function. This is a topic for future 
research, which requires a relatively sophisticated treat- 
ment of (23). 

5.5. Propagation 

The ability of low mode internal waves to propagate 
long distances from their generation sites is loosely con- 
strained in the observational record and is a serious is- 
sue for internal wave energy balance studies. The most 
significant processes are tidal constituents and near- 
incrtial waves. These internal waves do not play a major 
role in this Regional study of the high wavenumber / 
high frequency spectral variability. They do, however. 



act as sources for the internal wave continuum, whether 
it be through extreme scale separated interactions (Sec- 
tion S.l.l.c) or via a PSI decay process of tidal energy 
into near-inertial waves that then is linked to the con- 
tinuum (Section 5.1.2.1). 

We anticipate that variability in the tidal and in- 
ertial peaks is linked to variability in the continuum. 
Fu [1981] documents variability of the inertial peak in 
the Western North Atlantic and obtains the following 
Regional classification scheme: Class 1 spectra have 
prominent inertial peaks. Such spectra occur in records 
obtained from the Mid-Atlantic Ridge in the vicinity 
of 28° N and occur in conjunction with relatively white 
frequency spectra. We anticipate that the inertial peaks 
here are related to a PSI decay process (Section 5. 1.2. a). 
Class 2 spectra have less prominent peaks and are found 
(a) in the upper ocean at depths less than 2000 m, (b) 
in the deep ocean above rough topography off the Mid- 
Atlantic Ridge axis and away from 28° N and (c) in the 
deep ocean under the Gulf Stream. We anticipate issues 
of internal lee-wave generation and wave-mean interac- 
tions (Section 5.2) to be significant. Finally, Class 3 
spectra have the smallest inertial peaks and are found 
above smooth topography. Fu [1981] argues that Class 
1 and Class 2 inertial peaks are associated with local 
processes and Class 3 inertial peaks are consistent with 
propagation effects as waves initialized with the GM76 
spectrum refiect from their turning latitudes. Such 
quantitative estimates depend upon the vertical modal 
distribution of energy, and the regional variability of m^ 
is not well constrained in that study. Fu [1981] 's study 
included the MODE Central mooring instrumented at 
[500 (class 2a), 1500(class 2a), 4000(class 3)] m depth, 
whose location was the site of moorings during IWEX 
and FASINEX, Section 3.3.2; and PolyMode IIIc moor- 
ing 82 [172(class 2a), 322(class 2a) 522 (class 2a), 2446 
(class 2b), 3946 (class 2b)], Section 3.3.3. The vertical 
profile data sets at these locations in our study are fo- 
cused upon depths less than 2000 m. Leaman [1976], 
however, comments upon differences between vertical 
profiles obtained over rough and smooth topography in 
the vicinity of MODE Central. 

The last decade has seen the development of Gen- 
eral Circulation Models that include both tidal forcing 
and wind products with sufficient temporal resolution 
to force inertial motions [Simmons, 2008; Arbic et al., 
2010]. The models, however, obey their own dynam- 
ics: near-inertial forcing will be influenced by the mixed 
layer scheme, absence of an eddy field, presence/absence 
of a realistic internal wavefield, topographic representa- 
tion, sub-grid closures, etc. Unraveling how such model 
fields relate to the ocean will be a complicated, but use- 
ful aid in interpreting the oceanic internal wavefield. 

The issue of propagation demands careful assessment. 
As long as the roles of the bottom boundary condition 
and refraction in the mesoscale eddy field are not ap- 
preciated, the dominant balance of terms in an energy 
budget can seriously be misconstrued. See, for exam- 
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pie, the difference of opinion about the Brazil Basin 
data set in St. Laurent and Garrett [2002] and Polzin 
[2004b] that appears to be resolved in favor of a high 
efficiency of topographic scattering in Polzin [2009] . 

Similarly, Afford and Zhao [2007] estimate significant 
horizontal energy ffuxes in the near-inertial band north 
of the Gulf Stream and regard those fiux estimates as 
symptomatic of the ability of low-mode waves to propa- 
gate large distances. The issue of mooring motion con- 
taminating the horizontal energy flux estimates bears 
further investigation before those flux estimates can be 
considered robust. 

6. CONCLUSION 

The ingredients that shape the deep ocean energy 
spectrum have been known for some time. Fig. 34 from 
McComas [1977] presents one version of how these in- 
gredients can be combined. In that recipe, the primary 
sources are at high frequency and a transfer of energy 
to higher vertical wavenumber is associated with the 
ID mechanism. Parametric spectral fits summarized in 
Fig. 37 arc highly suggestive that the ID mechanism 
plays an active role in determining the observed power 
laws. However, the details are not consistent with such 
an interpretation. 

We now believe that near-inertial and tidal sources 
dominate high frequency sources of internal wave en- 
ergy, and this posses a real challenge to the induced dif- 
fusion mechanism: Energy transports associated with 
induced diffusion are to higher vertical wavenumber and 
lower frequency for vertical wavenumber domain power 
laws 6 > 0. That requires a source of high frequency 
energy. Possible sources are coupling of internal waves 
with surface swell, internal wave coupling to mixed layer 
turbulence and a weaker version of internal wave cou- 
phng to the eddy field. These possible sources are not 
well constrained, either observationally or theoretically. 

A second facet is that transports under the resonant 
interaction approximation are not completely under- 
stood. Hitherto neglected local interactions may play 
a significant role in determining nonlinear transports. 

A third facet is that resonant kinetic equation eval- 
uations in Fig. 40 indicate GM76 is far from being a 
stationary state, seemingly at odds with GM76 repre- 
senting an ID no flux state. Self-consistent solutions to 
a broadened kinetic equation could be much closer to 
being stationary and exhibit energy transfers appropri- 
ate for low frequency sources. A closely related result 
is that ray-tracing diagnostics imply a transport of ac- 
tion to higher frequency. This transport could provide 
a pathway to supply energy to the high wavenumber 
continuum. A flrst principles derivation of the spec- 
tral transports associated with the ray-tracing models 
would be enlightening. 

Despite such uncertainty, nonlinearity is clearly an 
organizing principle. Power laws associated with para- 



metric spectral representations lie close to the Induced 
Diffusion stationary states of the resonant kinetic equa- 
tion describing the lowest order nonlinear transfers. 
The one exception appears to be a wavefleld set up by 
the decay of a semidiurnal internal tide through the 
parametric subharmonic instability. 

While our ideas of the ingredients and their geo- 
graphic distribution has evolved, our knowledge of how 
to create the observed spectrum (the recipe) has re- 
mained relatively static. We are confldent that such 
recipes are within reach, but it will take a combina- 
tion of observations, analytic work in the context of a 
radiation balance equation and realistic direct numeri- 
cal simulations. Theoretical studies coupled to observa- 
tional programs, especially those that would document 
2-D spectra, transfer rates, the spectral character of 
the forcing functions and seek to define the role of the 
geostrophic flow fleld are a high priority. Resonant cal- 
culations show strong non-stationary tendencies, and 
this issue needs to be either resolved or circumvented 
before addressing the contributions of other ingredients 
to the recipes. 

Regionality 

Variability in the spectral characteristics of the deep 
ocean internal wavefleld documented in Section 3 and 
consideration of regional variability in the forcing fields 
in Section 5 leads us to propose the following strawmen 
for further investigation: 

• The Natre region is a minimum for eddy energy in 
the North Atlantic and wind work at near-inertial fre- 
quencies is small. Equatorward of the critical latitude 
for M2 tides, PSI decay of the low-mode internal tide 
represents the major forcing of the internal wavefield. 
The Natre spectrum is an outlier with extremely high 
bandwidth, steep vertical wavenumber and shallow fre- 
quency spectral slopes. A preliminary recipe for the 
Natre region is described in Section 5. 1.2. a. This recipe 
invokes a heuristic cascade closure for nonlinearity that 
has its roots in the ray-tracing diagnostics. 

• The contrasts between the waveflelds north and 
south of the Gulf Stream are distinct and appear to have 
a parallel with differences between subtropical and sub- 
polar gyres in the North Paciflc along its eastern bound- 
ary. The later case deserves further investigation to 
confirm spectral parameters graphically extracted from 
the literature. Our leading hypothesis is that the differ- 
ences can be directly related to the relative roles of near- 
inertial wind forcing and interactions with the mesocale 
eddy fleld via (71) and (72). We have not yet translated 
this hypothesis into a recipe. 

• The tentative signature of a phase lag between 
seasonal cycles in high frequency and near-inertial fre- 
quency waves in the eddy desert of the North Paciflc 
may be a signature of nonlinear transfers in the absence 
of wave-mean coupling. 

The future will certainly bring many further exciting 
developments, and the synthesis of theoretical, observa- 
tional and numerical results yet to be obtained. 
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TABLE 2. Meta data for moored current meters used in this study. 



place 


latitude 


longitude 


water depth 


sensor depth 


inat. type 


code 


start 


atop 


N cph 


0, C/m 


Sitc-D 


39 17.9' 


70 05.6' 


2640 m 


106 
511 
1013 
1950 


VACM 


WHOI 2203 
WHOI 2204 
WHOI 2205 
WHOI 2206 


Feb. 26, 1967 


Mar. 25, 1967 


3.03 
1.48* 
0.66* 
0.58* 




Shelf- Slope 


39 36.6' 


70 56.5' 


2305 m 


305 


VACM 


WHOI 5881 


Feb. 10, 1976 


Aug. 8. 1976 


2.82* 


-0.0235* 


Primer 


39 05.2' 


69 21.4' 


2990 m 


391 


VACM 


WHOI 9872 


Dec. 8, 1995 


Dec. 4. 1997 


2.39* 


-0.0170* 


Iwex 


27 43.9' 


69 50.95' 


5453 m 


604 


VACM 


WHOI 515A4 


Nov. 3, 1973 


Dec. 15, 1973 


2.55 


-0.0175 










604 


VACM 


WHOI 515B4 


Nov. 3, 1973 


Dee. 15, 1973 


2.55 


-0.0177 










633 


VACM 


WHOI 515C5 


Nov. 3, 1973 


Dee. 15, 1973 


2.60 


-0.0157 


Lotus 


34 1.2 


70 1.45 


5366 


328 


VACM 


WHOI 7666 


May 11, 1982 


April 12, 1983 


1.49 


-0.00356 




33 58.55 


70 0.36 


5366 


348 


VACM 


WHOI 7886 


April 14, 1983 


May 1, 1984 


1.49 


-0.00356 


PolyMode III-c 


15 23.40 


53 55.20 




319 


VACM 


NOVA 80 


May 11, 1977 


May 1, 1978 


2.78 


-0.03096 




15 11.50 


53 12.30 




309 


VACM 


NOVA 81 


May 12, 1977 


Dee. 20, 1977 








15 02.10 


54 12.90 




338 


VACM 


NOVA 82 


May 13, 1977 


May 1, 1978 






Fasincx 


27 58.9 


69 58.8 




556 
631 


VACM 
VACM 


WHOI 8293 
WHOI S294 


Oct. 29, 1984 


June 18, 1986 


2.48 
2.61 


-0.0178 
-0.0202 


Subduetion 


25 31.90 


28 57.2 




300 


VMCM 


WHOI 


Feb. 12, 1992 


Oct. 14, 1992 


2.44 


-0.02244 


Subduetion 


25 31.90 


28 57.2 




300 


VMCM 


WHOI 


Oct. 15, 1992 


•June 16, 1993 


2.44 


-0.02244 


EBC 


37 6.7 


127 32.1 


4752 


598 


Aandcraa 


rcm02674 


Aug. 9. 1992 


Aug. 19, 1994 


2.01 


0.00455 


Ocean Storms 


47 25.4 


139 17.8 


4224 


199 


VMCM 


rcm07420 


Aug. 20, 1987 


Juno 25, 1988 


2.25 


0.006 


AIWEX 


74 


143-144 


3700 


250 


SBE-3 




Mar. 20, 1985 


Apr. 5, 1985 




0.02 



APPENDIX A: META DATA 

Meta data for the data sets used in this study can be 
found in Tables 2 and 3. 



APPENDIX B: INSTRUMENTATION 

A terse presentation of many of these data was made 
in Lvov et al. [2004], simply summarizing the avail- 
able power law estimates as points in a frequency- 
wavenumber domain. That presentation was limited 
to spectra appearing in published literature. The in- 
tent of this work was to present pertinent data in a 
common framework, limiting the potential pitfalls asso- 
ciated with an irregular analysis. This does not limit 
uncertainties associated with the use of different instru- 
ment systems. Details of an instrumental and technical 
nature that impact the interpretation of the data are 
collected herein. 

Bl. Moored Current Meters 

Regarding the interpretation of moored current me- 
ter data as internal waves, there are two primary sets 
of issues. The first set is referred to as finestructure 
contamination. Here the presumption is that the data 
record represents signal and that departures from lin- 
ear internal wave kinematics are associated with either: 
(a) quasi-permanent finestructure (either density or ve- 
locity) being advected past the sensor (e.g., Polzin et 
al. [2003], 2.1.4), (b) self-advection within the wavefield 
(e.g., Sherman and Pinkel [1991], 5.1.4), or (c) attempt- 
ing to estimate density perturbations using temperature 
data only and invoking a stable relation between den- 
sity, temperature and salinity. This later assumption 
fails in regions that have significant large-scale gradients 
of temperature and salinity on isopycnals (e.g., Ferrari 
and Polzin [2005]). See the cited references for further 
discussion of these issues. The second set of issues are 
instrumental in nature: (a) drag associated with flow 
past a mooring will induce movement (e.g., FofFonoff 
[1966]), (b) typical sampling rates are not sufficient to 
resolve internal wave spectrum out to the buoyancy fre- 



quency in the upper ocean, and (c) instrument response 
issues. 

Regarding (a) , there are typically 3 classes of moor- 
ings: deep, intermediate and surface. A deep mooring 
will utilize glass balls as buoyancy elements distributed 
along the mooring cable. The primary buoyancy ele- 
ment for an intermediate mooring will be a large syn- 
tactic sphere at the upper terminus. This mooring type 
utilized at WHOI starting in the early 1980's. It has less 
over all drag and consequently provides a more stable 
platform. Finally, surface moorings are loosely tethered 
in order to accommodate the often sizable surface mo- 
tions associated with surface waves. As a result, the 
surface buoy inscribes a watch circle roughly equal to 
the water depth [Trask et al., 1982]. Vertical motion 
of the instruments can be diagnosed if pressure sensors 
are included, e.g., Fig. 51. Experience suggests iner- 
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Figure 51. LOTUS frequency spectra of horizontal ki- 
netic energy (blue line), potential energy (green line), 
and pressure p interpreted as a vertical displacement, 
N^p^ /2. This interpretation implies possible contamina- 
tion by mooring motion is more problematic at inertial 
and tidal frequencies than within the internal wave con- 
tinuum. 
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TABLE 3. Meta data for the vertical profiling instrumentation used in this study 



place 


latitude 


longitude 


water depth 


inst. type 


start 


stop 


Station-W{3) 


39 


70 


3000 


MMP 


Jun. 20, 2003 


Aug. 5, 2003 


Station-W(3) 


39 


70 


3000 


MMP 


Jan. 17, 2004 


Mar. 2, 2004 


Mode 


28 


69 40' 


5440 


EMVP 


Jun. 11, 1973 


Jun. 15, 1973 


Fasinex 


27-29 


67 30' - 70 30' 


S 5200 


HRP 


Feb. 17, 1986 


Mar. 5 1986 


Natre 


23 56' - 27 32' 


26 45' - 30 43' 


S 5500 


HRP 


Mar. 28, 1992 


Apr. 14, 1992 


SFtre2 


14-16 


50-57 


S 5400 


HRP 


Nov. 14, 2001 


Nov. 24, 2001 


Patchexl 


34 


127 


S 4700 


MSP 


Oct. 17, 1986 


Oct. 24, 1986 


Patchex^ 


34 


127 


S 4700 


FLIP 


Oct. 1986 


7.5 days 


Swapp 


35 8.2' 


126 59.0' 


S 4700 


FLIP 


Feb.-Mar. 1990 


18 days 


STREX 


50 


140 


S 4000 


XCP 


Nov. 8 1980 


Nov. 22, 1980 


AIWEX 


74 0' - 74 12' 


144 0' - 145 0' 


S 3750 


XCP 


April 3, 1985 


April 14, 1985 



tial and tidal frequencies suffer greater contamination 
by mooring motion than super-tidal frequencies. Signif- 
icant discussions appear in Fu [1981] and Ruddick and 
Joyce [1979]. 

Regarding (b), insufficient sampling rates rates will 
result in the aliasing of high frequency signals back 
into the resolved frequency domain. Ordinarily, the fre- 
quency spectra will be sufficiently red that the aliased 
energy is relatively small. However, if there is a sub- 
stantial bump at the buoyancy frequency associated 
with turning point dynamics (e.g., Desaubies [1975]), 
the aliased energy may be significant. Ascertaining 
whether this is the case requires information on the ver- 
tical wavenumber content of the wavefield and details of 
the vertical structure of the buoyancy profile. Without 
such information and in the absence of pressure records, 
it is difficult to assert that the departure of the observed 
high frequency spectra in Fig. 5, for example, are noise 
rather than signal. 

Regarding (c), the principle moored current me- 
ters used here are the Vector Averaging Current Me- 
ter (VACM) and Vector Measuring Current Meter 
(VMCM). The VACM employs a Savonius rotor and 
the VMCM uses sets of propellors. The VMCM [Weller 
and Davis, 1980] was designed to have a cosine response 
to eliminate 'pumping' associated with the motion of a 
surface mooring. VACMs are standard for subsurface 
moorings. VACMs have a finite stall speed of about 2 
cm s~^ [Luyten and Stommel, 1991]. Good directional 
data is believed to be obtained at speeds smaller than 
this. The resulting 'noise' is not well defined. VMCM 
data are calibrated assuming no significant stalling. See 
Hogg and Frye [2007] for a more recent discussion of 
current meter performance. 

B2. Doppler Sonars 

Doppler sonars do not measure horizontal velocity as 
a point but determine Doppler shifts as a function of 
range along multiple acoustic beams. These Doppler 
shifts are interpreted as slant velocities of suspended 
acoustic backscattering targets moving with the water 
relative to the instrument platform. Individual esti- 
mates of ocean velocity are highly uncertain: hundreds 
of pings are typically used to average incoherent noise 
from the system. Estimates of signal-to-noise ratios ap- 



pear in Anderson [1992] and an appendix of Alford and 
Pinkel [2000]. The insidious complication of package 
motion and beam separation effects are considered in 
Polzin et al. [2002]. 

B3. Vertical Profiling Instrumentation 

Most of the vertical wavenumber domain information 
considered here was obtained with vertical profiling in- 
strumentation. Estimates of the vertical structure of 
the velocity field can be obtained by using a number of 
different sensors, each of which with its own strengths 
and weaknesses. 

Sanford's ElectroMagnetic Velocity Profiler (EMVP, 
Sanford [1975]) senses the voltage drop associated with 
an electrical current (seawater) in a magnetic field (the 
Earth's, in this case). The measurement is uncertain 
to within a conductivity weighted mean, so that other 
sensors need to be incorporated to provide estimates 
of the absolute velocity field. Noise levels are typically 
around 0.5 cm s~^, which is sufficient to start resolv- 
ing vertical wavelengths smaller than 10m. Expendable 
Current Profilers (XCPs) operate on the same princi- 
ple, but may have somewhat larger noise levels (1.0 cm 
s-i). 

Data from several acoustically tracked dropsonde's 
have been analyzed but not presented. With this 
instrument system, the horizontal position is esti- 
mated from range information provide by an acous- 
tical net and depth information from an onboard 
CTD (Conductivity- Temperature-Depth sensor suite). 
Oceanic velocity estimates assume that the package 
tracks the flow field as it descends. This method pro- 
vides estimates of the absolute velocity, but the small 
scale noise is relatively high. Acoustic tracking provides 
positions to within 1 m. A typical descent rate of 1 m 
s~^ implies that the background internal wavefield hav- 
ing vertical wavelengths smaller than several hundred 
meters will be obscured by noise. 

A third method is to utilize an acoustic travel time 
sensor to estimate relative flow past a freely-falling ve- 
hicle [e.g., the High resolution Profiler [Schmitt et al., 
1988] and the Multi-Scale Profiler [Winkel et al. , 
1996]]. Estimates of the oceanic velocity profile are de- 
duced from these relative velocity data and a model 
of how the vehicle responds to the relative flow. The 
method is capable of resolving oceanic shear at 1-m 
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scales. The limitations are at larger wavelengths. The 
offsets (zeros) of the acoustic travel time sensor are typ- 
ically determined in situ and may even be temperature 
and pressure dependent. Thus the resulting profiles are 
uncertain to a linear trend. Other information, such as 
provided by an electric field sensor or matching with 
a shipboard ADCP record, are required to provide the 
largest scale information. 

Acoustic travel time sensors are used in the McClane 
Moored Profiler (MP, Doherty et al. [1999]). This in- 
strument autonomously samples relative velocity, tem- 
perature, conductivity and pressure while transiting a 
mooring cable. Estimates of relative velocity are thus 
absolute. Noise levels in the velocity record are nomi- 
nally estimated at 0.5 cm s"^. Biases related to drifting 
zeros and mooring oscillations are not presently well de- 
fined. 

Velocity data obtained with Doppler sonars typically 
estimate horizontal velocity from back-to-back acoustic 
beams. In so doing, one assumes that the velocity field 
is horizontally uniform. This assumption breaks down 
as the beam separation increases and will preferentially 
contaminate high frequency, small vertical scale signals. 
Polzin et al. [2002]. 

APPENDIX C: PROCESSING 

In the absence of a mean flow, the linear internal 
wave equation is (e.g.. Gill [1982]): 

(9,2 + f)dlw + [N%z) + d^,]{dl + dl)w ^ (CI) 

for arbitrary stratification profile N'^{z) and vertical ve- 
locity w. If the stratification profile varies much more 
slowly than the wave phase, a WKB approximation for 
vertically propagating waves provides the approximate 
solution: 



w oc 



iV(z)-i/2e'/^(^)'^^ 



(C2) 



and so the effects of a variable buoyancy profile can be 
accounted for by stretching the depth coordinate by N 
and scaling the horizontal velocities by N^^^/No , in 
which No is a reference stratification. The value No — 
3 cph is often used. The use of the WKB approxima- 
tion requires w^ <C iV^. If this relation is not satisfied, 
solutions can be found by treating (CI) as an eigen 
value problem with appropriate boundary conditions. 
If cj^ <C iV^ but the wave-phase is not slowly varying, 
the boundary conditions arc that w — ai the top 
{z = 0) and bottom (z = H), which then implies the 
horizontal velocities [u{z),v{z)] are proportional to: 



N{zy'^cos{mT / N{z')dz'/ / N{z')dz'), (C3) 



H 



for integer values of n. For data sets that document 
the velocity profile over the entire water column, the 
first three modes have been estimated by using a linear 



regression. Vertical wavcnumber spectra are calculated 
by using a cubic spline to interpolate onto a uniform 
grid having approximately the same resolution as the 
unstretched data. 

Current meter data used here have sampling inter- 
vals of 3.75 minutes to 1 hour. Frequency spectra have 
been calculated using transform intervals of 10| and | 
days. Spectra are displayed using the later at super- 
tidal frequencies {a > 2.4 / cpd) and the former at 
lower frequencies. The intent is to give finer resolution 
to inertial and tidal peaks and decreased variability at 
higher frequency. 



APPENDIX D: A SARGASSO SEA SPECTRUM 

We have characterized the background internal wave 
spectrum in the Sargasso Sea as: 



E{m, a) 
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with ?Ti* ~ 107r/1300 m^^ and normalization con- 
stant A chosen to return E — 0.56 x 30 x lO"** m^ 
s~2 . This places the gradient variance at approximately 

3 

4 



p27r/10n 



L"' '^"'"^ 2m^E{m)dm — '^N'^- Power laws are char- 
acteristic of the Sargasso Sea Data sets (Section 3.3.2) 
are selected specifically to lie on the Induced Diffusion 
constant fiux solution (Section 5.1.1.c). The inertial 
peak was selected with the constraint that Ek/Ep = 4.4 
based upon moored data. We have deferred from invok- 
ing non-separable conditions via m^ — m^(a) that are 
a documented part of the observational record [Miiller 
et al., 1978] and appear as Ek{m — 0.01cpm)/i<^p(m = 
G.Olcpm) in excess of this value [Polzin et al., 2003]. 

ACKNOWLEDGMENTS. 

The authors gratefully acknowledge the efforts of many tal- 
ented researchers who were instrumental in collecting, pro- 
cessing and archiving the data presented here. Special thanks 
are extended to John Toole for providing yet unpublished MP 
data from Site-D, Tom Sanford for providing access to a num- 
ber of his EMVP and AVP data sets, and to Carl Wunsch 
for the PolyMode IIIc current meter data set. We greatfuUy 
acknowledge funding provided by a Collaborations in Mathe- 
matical Geosciences (CMG) grant from the National Science 
Foundation. 



NOTES 



1. Estimates of vertical bandwith jt can be obtained by assuming 
a parametric spectral representation and fitting lagged coher- 
ence estimates to the inverse Fourier transform of the spec- 
trum, [Garrett and Munk, 1972; Cairns, 1975; Miiller et al., 
1978] . This technique can be found in many textbooks as the 
Wicner-Khinchine relation. 
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2. Departures from a simple power law behavior arc apparent 
in the thermocline data at high frequencies and spectra from 
the deeper current meters do not roll-off at frequencies exceed- 
ing the local TV. Data from the very stable trimoored IWEX 
mooring roll-off nicely (see below) and thus Briscoe [1975b] 
suggests the super-buoyancy extension is an artifact of moor- 
ing motion. On the other hand, this does not mean that the 
departures from a simple power law behavior apparent in the 
upper current meters is an artifact of mooring motion. See 
Appendix Bl for further details. 

3. One also finds mud waves. Mud waves are sedimentary features 
of 1-10 km horizontal wavelength having amplitudes of lO's to 
lOO's of meters. These horizontal scales arc appropriate for the 
generation of freely propagating internal lee waves (with Eule- 
rian frequency cr = 0) if the intrinsic frequency cr = cr— p • u lies 
between the Coriolis and buoyancy frequencies: / < p ■ u < A'^. 
Significant coupling between the 'mean' and internal wavefield 
is anticipated at mean flow rates of 0.1-0.2 m s~^ . Sediment 
transport is an issue at such flow rates and the possibility ex- 
ists that the lee wave velocity perturbations affect the deposi- 
tion/erosion process so as to reinforce the mud-waves [Blum- 
sack, 1993]. But this gets us off the topic of the background 
internal wavefleld. 

4. The scale invariance assumption is, formally, a generic assump- 
tion about dependent and independent variables in the kinetic 
equation. Sec, for example. Bender and Orzag [1978], their 
Section 1.7 

5. A contrary assessment is given in Lvov et al. [2004]. That re- 
sult is associated with a sign error in the numerical evaluation 
of (37). 

6. Hendry's estimate of 300 W/m to the south east represents a 
short term estimate from the MODE current meter array 
centered about 28°N, 70°W. Long term estimates from sin- 
gle moorings at 31°N, 70°W (LDE) and 34°N, 70°W (LO- 
TUS) [Alford and Zhao, 2007] are approximately 150 W/m 
to the north-northwest and 200 W/m to the southeast. It 
is tempting to interpret energy flux estimates in Fig. 1 of 
Simmons [2008] as suggesting that the western boundary is 
a source of baroclinic tidal energy at these latitudes. The 
lack of consistency is troubling and the issue warrants further 
investigation. 
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